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ABSTRACT 


Studies  of  the  electron-density,  electron  and  ion  temperatures  in  the  F-region  were  made  by  means 
of  ground-based  radar  observations  at  the  Millstone  Hill  Radar  Observatory.  A  70-meter  parabolic 
antenna  directed  vertically  and  a  2.  5-Mw  pulse  radar  operating  at  440  Mcps  were  employed  for 
these  measurements  which  were  conducted  for  periods  of  30  hours  at  approximately  weekly  intervals 
throughout  1963.  Examination  of  the  echo  power  as  a  function  of  height  leads  to  a  profile  of  elec¬ 
tron  density  with  height,  provided  the  electron  and  ion  temperatures  are  the  same  (T  =  T.).  Addi¬ 
tional  measurements  of  the  spectra  of  the  signals  corresponding  to  different  heights  permit  the  ratio 
T^/T.  to  be  determined  and,  where  this  is  different  from  unity,  the  observed  profile  can  then  be 
corrected  for  the  effect  on  the  scattering  introduced  by  the  inequality  in  temperature. 

Results  of  observations  extending  over  a  period  of  one  year  from  February  1963  to  January  1964  are 
presented.  The  ratio  T  /T.  achieved  a  maximum  value  ~2. 0  to  2.  6  at  a  height  of  about  300  km 
soon  after  dawn,  irrespective  of  the  season.  There  was  little  change  in  height  dependence  in  this 
ratio  throughout  the  daylight  hours,  and  at  sunset  the  ratio  fell  with  a  time  constant  of  the  order  of 
an  hour.  At  night  T  /T.  was  occasionally  close  to  unity,  but  more  often  a  significant  difference 
remained  in  the  temperatures  at  all  heights.  Ion  temperature  increased  with  height  at  all  times, 
but  above  500  km  this  may  be  due  in  part  to  the  presence  of  an  unknown  amount  of  He+  ions,  which 
considerably  affects  the  interpretation  of  the  signal  spectra.  Electron  temperatures  were  largely 
independent  of  height  above  about  300  km.  Evidence  is  presented  of  ionospheric  heating  during 
magnetically  disturbed  conditions,  but  it  is  shown  that  this  is  only  of  great  importance  at  night. 
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IONOSPHERIC  BACKSCATTER  OBSERVATIONS 
AT  MILLSTONE  HILL 


I.  INTRODUCTION 

In  recent  years  the  study  of  the  upper  part  of  the  earth’s  ionosphere  has  received  considerable 
impetus  from  the  development  of  satellite  and  rocket  methods  of  exploration.  In  the  measurement 
of  electron  density,  outstanding  results  have  been  obtained  with  the  satellite  Ariel  and  the  top¬ 
side  sounder  Alouette.  In  addition,  the  rocket  measurements  made  from  Wallops  Island,  Vir- 

1 

ginia  have  been  particularly  successful.  Partly  as  a  result  of  this  work,  and  partly  from  the 

study  of  the  drag  on  satellites  produced  by  the  neutral  atmospheric  constituents,  early  ideas 

concerning  the  behavior  of  this  region  have  been  abandoned.  Until  very  recently  it  was  widely 

believed  that,  above  the  peak  of  the  F-region,  0+  was  the  principal  ion  over  a  wide  range  of 

heights.  It  was  also  thought  that  the  temperatures  of  the  ions  T\  and  electrons  were  equal 

and  independent  of  height,  and  that  the  electron-density  distribution  was  wholly  controlled  by 
2-5  • 

diffusion.  Thus  the  electron  density  would  decrease  with  altitude  with  a  scale  height 

H  =  k(T  +  T.)/m.g  =  2kT./m.g,  where  k  is  Boltzmann’s  constant  and  m.  is  the  mass  of  atomic 

oxygen.  In  this  model,  H  is  not  expected  to  vary  significantly  with  height,  and  some  rocket 

experiments  supported  this  prediction.  This  model  has  given  way  to  one  more  complex  as  a 

consequence  of  observations  which  show  that:  (a)  the  uppermost  parts  of  the  ionosphere  con- 

6-9 

sist  of  regions  in  which  helium  and  hydrogen  ions  predominate,  (b)  the  electron  and  ion 

10-14  +  + 

temperatures  need  not  be  the  same,  (c)  charge  exchange  between  O  and  H  (Ref.  15) 

governs  the  proton  distribution,  and  diffusive  equilibrium  of  the  separate  constituents  may  not 
be  established  until  heights  of  700  km  are  reached.^ 

The  results  obtained  by  the  topside  sounder  Alouette  indicate  that  very  rarely  does  the 
upper  F-region  electron  density  exhibit  constant  scale  height  H  (Refs.  17  and  18).  This  result 
is  interpreted  as  indicating  that  one  or  more  of  the  above  effects  operate  at  most  latitudes.  In 
addition,  this  satellite  has  shown  that  the  equatorial  anomaly  observed  with  ground-based  iono- 
sondes  extends  into  the  upper  F-region  over  a  range  of  geomagnetic  latitudes  of  at  least  ±20° 
(Refs.  19  and  20).  Over  these  regions,  and  possibly  others,  the  geomagnetic  field  plays  a  signif¬ 
icant  role  in  governing  the  electron-density  distribution.  Theoretical  attempts  to  account  for 

21-24 

this  control  have  recently  been  made,  but  it  cannot  yet  be  claimed  that  all  the  features  are 

19 

understood. 

At  high  latitudes  geomagnetic  control  may  disappear  but  significant  fluxes  of  precipitating 
particles  may  contribute  to  the  ionization,  particularly  in  the  auroral  zones.  It  seems,  there¬ 
fore,  that  only  in  temperate  latitudes  is  the  F-region  likely  to  be  largely  uninfluenced  either  by 
geomagnetic  control  of  diffusion  or  additional  sources  of  ionization  from  the  Van  Allen  belts. 
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In  these  latitudes  an  understanding  of  the  behavior  of  the  ionospheric  F-layer  could  perhaps  be 
derived  from  observations  in  which  the  following  quantities  are  measured  simultaneously  as 
functions  of  height:  (a)  electron  density;  (b)  electron  temperature;  (c)  ion  temperature;  (d)  ion 
composition;  (e)  neutral  density;  (f)  neutral  temperature;  (g)  neutral  composition;  and  (h)  the 
incident  solar  radiation.  At  the  present  time  many  of  these  quantities  are  being  measured  in¬ 
dividually  by  means  of  satellites,  but  it  is  often  difficult  from  these  data  to  extract  the  separate 

dependence  of,  for  example,  the  electron  temperature  on  height,  local  time,  latitude  (and  longi- 

25 

tude)  and  variations  in  the  solar  flux.  Thus,  rockets  appear  to  provide  the  best  means  of  mak¬ 
ing  these  measurements,  but  so  far  none  have  been  equipped  to  measure  simultaneously  more 
than  three  of  the  quantities  involved.  Further,  in  order  to  observe  temporal  variations,  repeated 
rocket  launchings  are  required  involving  considerable  effort  and  expenditure. 

One  relatively  inexpensive  ground-based  method  can  provide  measurements  of  the  quantities 

(a)  through  (d),  but  being  a  radio  method  it  is  insensitive  to  the  quantities  (e)  through  (h).  This 

26 

method  is  the  radar  incoherent  backscatter  technique  discussed  by  Gordon  and  employed  to 

27  28  29  30 

study  the  ionosphere,  first  by  Bowles  '  and  later  by  Pineo,  et  al.  '  The  backscatter 
facility  in  its  present  form  came  into  operation  in  January  1963,  and  this  report  presents  the 
results  of  the  first  year  of  operation.  In  Sec.  II  we  review  briefly  the  theory  of  the  method. 

A  detailed  account  of  the  incoherent  backscatter  radar  facility  at  the  Millstone  Hill  Radar 
Observatory  is  provided  in  Sec.  Ill,  the  observational  methods  are  presented  in  Sec.  IV,  and 
reduction  procedures  are  presented  in  Sec.V. 

II.  THEORY 

A.  Electron -Density  Distribution 

26 

In  1958,  Gordon  pointed  out  that  if  a  very  powerful  radar  beam  were  directed  at  the 
ionosphere,  a  weak  but  detectable  echo  from  the  free  electrons  there  might  be  obtained.  The 
intensity  of  the  echo  can  be  computed  by  assuming  that  the  electrons  scatter  independently  with 
a  scattering  cross  section  ac  given  by  the  square  of  the  classical  radius  rg  =  (e2/mc2).  Some 
confusion  has  arisen  because  this  cross  section  is  defined  as  that  which  scatters  energy  into 

unit  solid  angle,  whereas,  in  radar  calculations,  it  is  customary  to  normalize  the  cross  section 

31 

to  correspond  to  power  reflected  into  47 r  solid  angle.  Hence  the  radar  cross  section  a  is 

-28  2  ^ 
expected  to  be  47ro-c  (~10  m  ).  For  a  volume  containing  N  electrons,  the  phase  of  the  N  re¬ 
flected  waves  will  be  independent  and  the  powers  should  add  to  give  an  average  scattered  power 

proportional  to  Ncr  .  Because  at  any  given  height  the  electrons  completely  fill  the  radar  beam, 

e  /  2  /  4 

the  echo  power  should  vary  with  range  R  only  as  l/R  ,  not  1/R  as  for  conventional  (i.e., 

” point”)  targets.  For  observations  conducted  in  the  zenith  it  may  be  shown  that  the  echo  power 
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is  given  in 


(1) 


P 


r 


where  P^  is  the  transmitter  power,  tj  ^  is  the  efficiency  of  the  waveguide  or  feeder  system,  c 
is  the  velocity  of  light,  r  is  the  pulse  length,  X  is  the  radio  wavelength  and  h  is  the  height  to 
the  pulse  center.  G(0)  is  the  antenna  gain  over  a  lossless  isotropic  radiator,  assuming  that  the 
antenna  pattern  is  cylindrically  symmetrical  about  its  axis,  and  0  is  the  angle  subtended  in  any 


2 


direction  from  this  axis.  The  electron  cross  section  a  is  usually  different  from  or  =  (4ttc t  )  for 

33  e  c 

reasons  which  will  be  given.  For  a  typical  parabolic  antenna,  Evans  has  shown  that  Eq.  (1) 
can  be  reduced  to 


0.74 


P.r?  crNaA 
t  'r  o 
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where  Aq  is  the  effective  antenna  aperture.  From  Eq.  (2)  we  see  that,  for  the  detection  of  a 
given  density  N  at  a  given  height  h,  the  radar  designer  can  adjust  only  the  product  P^tA^.  If, 
however,  it  is  necessary  to  obtain  a  given  height  resolution  Ah  =  ct/2,  one  has  no  choice  but  to 
employ  both  a  high  transmitter  power  P^  and  large  effective  aperture  Aq  to  obtain  useful  results 
One  of  the  largest  systems  so  far  constructed  (at  Arecibo,  Puerto  Rico)  employs  the  same 

transmitter  power  P.  as  the  Millstone  Hill  radar,  but  has  a  value  of  A  which  is  20  times  that 

35  1  ° 

at  Millstone.  If  the  electron  density  were  independent  of  height,  the  Arecibo  radar  could  be 

expected  to  detect  echoes  four  times  higher  than  Millstone,  but  in  practice  this  ratio  will  be 

more  nearly  two  times. 


B.  Signal  Spectra 

28 

Bowles  was  first  to  discover  that  the  Doppler  width  of  the  reflected  signals  was  consider¬ 
ably  less  than  expected  for  a  model  of  the  scattering  in  which  it  is  assumed  that  each  electron 
contributes  a  signal  with  a  Doppler  shift  proportional  to  the  component  of  its  velocity  in  a 

direction  to  the  radar.  Bowles  correctly  attributed  this  behavior  to  the  role  played  by  the  ions 

3  6>  37—39 

in  the  plasma.  Subsequent  theoretical  papers  by  Fejer  and  others  7  have  shown  that  when 

the  radar  wavelength  X  is  greater  than  47rX^,  where  X^  is  the  Debye  length  (=  vkT/47j*Ne  ),  the 

simple  model  above  must  be  replaced  by  one  in  which  the  collective  behavior  of  the  electrons 

and  ions  is  considered.  The  motions  of  the  ions  govern  fluctuations  in  the  electron  density 

(with  a  scale  of  X^),  and  these  fluctuations  constitute  small  changes  in  the  refractive  index. 

When  the  radar  wavelength  is  significantly  larger  than  the  size  of  these  fluctuations,  they  must 

be  thought  of  as  the  scattering  centers  and  not  the  individual  electrons.  The  situation  becomes 

more  complex  when  observations  are  conducted  in  directions  normal  or  nearly  normal  to  the 

magnetic-field  lines,  for  then  the  gyrorotation  of  the  ions  causes  the  ion  motions  to  be  ordered. 

40  41-44 

This  case  has  been  considered  by  Fejer  and  others,  and  experimental  verification  of  the 

4  5 

narrowing  of  the  spectrum  has  been  reported  by  Pineo,  et  al.  As  magnetic  effects  play  no 

part  in  observations  at  the  Millstone  Hill  observatory,  they  will  not  be  discussed  further  here. 

3  6 

In  the  theory  developed  by  Fejer  it  is  assumed  that  both  electrons  and  ions  have  Boltzmann 
energy  distributions  and  that  collisions  are  infrequent  and  may  be  neglected.  Then,  if  it  is 
further  assumed  that  only  singly  charged  ions  of  mass  rm  are  present,  the  shape  of  the  echo 
spectrum  takes  the  form  shown  in  Fig.  1.  It  is  clear  from  this  figure  that  the  echo  power  fre¬ 
quency  spectrum  is  a  function  of  both  the  electron  temperature  T^  and  ion  temperature  T^; 
hence,  by  determining  the  shape  of  the  spectrum,  the  ratio  T  /T^  can  be  determined.  If,  in 
addition,  rm  is  known  then  T^  can  be  obtained  from  the  total  width  of  the  spectra.  The  first 

observations  of  the  spectrum  shapes  from  which  these  quantities  were  extracted  were  reported 
11  12 

by  Evans  and  Pineo  and  Hynek.  Subsequently,  further  measurements  have  been  reported  by 
Evans  and  Loewenthal,4^  and  Evans.4^4^ 
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C.  Electron  Cross  Section 


The  cross  section  of  the  electrons  cr  is  proportional  to  the  area  under  the  curve  in  Fig.  1 
and,  as  such,  depends  upon  the  ratio  Te/T„  This  dependence  has  been  considered  by  Bunneman 
and  others.50,  51  Bunneman  gives  for  c r  the  expression: 


Thus,  for  very  short  wavelengths,  47rX^  »  X,  <r  -*■  cr^.  In  the  ionosphere,  at  heights  presently 
being  studied  by  this  technique,  the  Debye  length  XD  is  of  the  order  of  a  few  millimeters,  and 
for  these  radars  X  ^  47rX^.  When  X  »  4ttXd,  Eq.  (3)  becomes: 

a  =  ,re{l/[l  +  (Te/T.)]}  ,  (4) 

and  thus  the  highest  cross  section  that  can  be  expected  (under  conditions  of  thermal  equilibrium) 
is  o-  =  0.5  a  (=  0.499  X  10"28m2). 

If  T^/T.,  is  a  function  of  altitude  h,  the  electron-density  profile  N(h)  can  only  be  determined 
from  Eq.  (2)  if  the  height  dependence  of  c j  is  known.  This  complication  was  not  immediately 

recognized  and  some  published  backscatter  profiles  of  electron  density  are  now  thought  to  be  in 

52  34 

error.  This  difficulty  has  caused  Bowles,  operating  at  50Mcps,  and  Greenhow,  et  al., 

operating  at  398Mcps,  to  adopt  a  method  of  deriving  an  electron-density  profile  from  measure¬ 
ments  of  the  Faraday  rotation  of  the  radio  waves  in  place  of  the  echo  power.  This  technique  is 
difficult  to  implement  at  frequencies  above  lOOMcps  where  the  amount  of  rotation  is  small. 

In  the  measurements  to  be  described,  the  echo  power  has  been  obtained  as  a  function  of 
height,  and  N(h)  obtained  from  Eq.  (2)  by  determining  the  variation  of  cr  with  height.  This  is 
accomplished  by  simultaneous  measurements  of  the  signal  spectra  corresponding  to  different 
heights  from  which  T^/T^  can  be  obtained.  The  height  dependences  of  Te/T.,  and  T.  are,  of 
course,  of  interest  in  their  own  right,  and  much  of  the  report  is  devoted  to  observations  of 
these  quantities. 

D.  Mixtures  of  Ions 

As  the  ion  mass  m.  enters  the  equations  in  the  form  (m.)1/2,  it  seems  that  0+  and  N+  (or 

+  +  +  1  1 

and  NO  )  could  be  considered  for  all  practical  purposes  as  identical,  in  that  it  would 

not  be  possible  experimentally  to  distinguish  one  from  the  other.  Below  about  200  km  it  is  ex- 

53 

pected  that  significant  concentrations  of  several  of  these  ions  exist,  and  Petit  has  considered 
the  experimental  accuracy  required  to  determine  the  ratio  of  the  heavier  to  the  lighter  ions. 
Observations  conducted  at  the  Millstone  Hill  observatory  during  1963  are  restricted  to  a  height 
range  of  approximately  200  to  750  km.  The  lower  limit  is  set  by  ground  clutter  echoes,  and  the 

upper  limit  by  the  strength  of  the  signals.  Over  much  of  this  height  range  0+  is  expected  to 

54  55  + 

predominate.  '  However,  above  about  750  km  by  day  and  perhaps  500  km  at  night,  H  and 

+  56  57 

He  ions  become  important  constituents.  '  Their  presence  can  modify  the  scattering  con¬ 
siderably  by  changing  the  shape  of  the  signal  spectra.  Figure  2  shows  the  influence  of  increas¬ 
ing  the  amount  of  He+  in  an  0+/He+  mixture  on  the  shape  of  the  spectra  for  the  case  T  /T.  =  2.0. 

+  e  1 

It  is  clear  that  the  presence  of  large  amounts  of  He  (e.g.,  50  percent)  can  be  recognized  readily. 
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However,  if  only  small  amounts  of  He+  are  present,  the  spectra  might  be  interpreted  as  indi- 

r  ft 

eating  an  erroneously  high  value  of  T.  and  low  value  of  T  /T..  Moorecroft;  has  considered 

1  e  1  +  +  + 

how  the  quantities  Tg,  T\  and  the  ratio  between  the  percentages  of  He  ,  H  and  O  ions  might 
be  extracted  from  spectrum  observations.  In  principle,  it  is  possible  to  perform  such  an 
analysis  provided  that:  (1)  the  measurements  can  be  made  to  a  very  high  degree  of  precision; 
and  (2)  it  is  possible  to  remove  some  ambiguities  in  the  interpretation  by  observations  conducted 
at  higher  and  lower  altitudes.  In  practice,  where  pulsed  radar  observations  are  made  one  must 
consider  also  the  distortion  of  the  theoretical  spectra  by  the  finite  width  of  the  transmitter 
pulses,  and  in  any  radar  system  there  will  be  some  further  smoothing  due  to  the  finite  width 
of  the  filters  employed  in  the  spectrum  analyzer.  Considerable  effort  has  been  made  to  compute 
these  effects  for  the  Millstone  Hill  Radar  Observatory  (Sec.  V). 

III.  EQUIPMENT 

A.  General 

The  Millstone  Hill  Radar  Observatory  is  located  in  Westford,  Massachusetts  (71.5°W, 
4Z.6°N).  The  parameters  of  the  ionospheric  radar  are  listed  in  Table  I.  The  antenna  (Fig.  3) 
is  a  70 -meter  parabola  directed  vertically.  The  feed  system  is  a  circular  horn  coupled  to  a 
turnstile  junction.  A  pair  of  opposite  ports  are  connected  into  remotely  adjustable  shorts 
allowing  any  mode  of  polarization  to  be  radiated.  In  practice,  circularly  polarized  waves  of 
one  sense  are  transmitted  and  the  opposite  sense  received.  The  transmitter  is  coupled  to  an¬ 
other  port  of  the  turnstile  junction  by  a  waveguide,  and  the  receiver  is  similarly  coupled  to  the 
fourth  port.  The  transmitter  and  receiver  frequencies  are  controlled  by  synthesizing  all  the 
needed  oscillator  frequencies  from  a  single  frequency  standard. 

B.  Receiver 

The  first  amplifier  in  the  receiver  is  a  Zenith  electron-beam  parametric  amplifier  which 
is  synchronously  pumped  (at  880Mcps).  This  amplifier  normally  introduces  additional  noise 
due  to  an  unwanted  side-band  or  image.  However,  by  employing  a  pump  frequency  which  is 
precisely  twice  the  radar  frequency  (i.e.,  synchronously  pumping),  the  image  can  be  removed. 


TABLE  1 

PARAMETERS  OF  IONOSPHERIC  RADAR 

Antenna 

70-meter  fixed  parabola 

Polarization 

Circular 

Effective  aperture 

1600  m2 

Frequency  f 

440  Mcps 

Transmitter  power 

2.5  Mw 

Transmitter  pulse  lengths  t 

0. 1-,  0.5-  and  1.0-msec  pulses  used 

Pulse  repetition  frequency 

50  cps 

Receiver  bandwidth  b 

25  keps 

System  temperature  T^ 

~200°K  (monitored  continuously) 

Post  detector  integration 

~20  db 

5 


Unfortunately,  the  phase  information  of  the  signal  is  then  destroyed;  that  is,  it  becomes  impos¬ 
sible  to  observe  at  the  output  of  the  receiver  if  the  signal  spectrum  is  asymmetric,  for  the 
amplifier  now  operates  both  in  the  conventional  way  and  as  a  side-band  exchanger.  The  statis¬ 
tical  fluctuations  of  the  noise  at  the  output  terminals  of  the  receiver  are  also  changed  so  that, 

after  n  independent  samples,  the  rms  temperature  fluctuation  is  ±*/2T  /\fn  instead  of  ±T  /\fn 

s  s 

as  in  the  case  of  a  conventional  amplifier  with  system  temperature  T  .  These  disadvantages 
are  outweighed  by  the  rapid  recovery  of  the  amplifier  following  the  transmitter  pulse,  its  com¬ 
plete  stability,  and  the  fact  that  it  cannot  be  destroyed  by  moderate  amounts  of  transmitter 
power. 

C.  Echo  Power  Integration 

The  electron-density  profiles  are  obtained  from  measurements  in  which  the  echo  power  as 
a  function  of  height  is  determined.  This  is  accomplished  by  placing  in  the  receiver  a  filter 
whose  bandwidth  (25kcps)  is  sufficient  to  accommodate  all  the  frequency  components  of  the  re¬ 
flected  signal.  The  output  of  this  filter  is  connected  to  a  square-law  detector  which  provides  a 
voltage  proportional  to  the  input  power.  This  voltage,  in  turn,  is  sampled  by  a  digital  voltmeter 
at  intervals  of  200  psec  and  assigned  one  of  256  possible  levels.  The  digital  numbers  represent¬ 
ing  these  voltages  are  continuously  summed  in  the  CG  24  computer  to  form  a  plot  of  the  integrated 
echo  power  vs  delay  over  the  time-base.  This  integration  process  is  continued  usually  for  periods 
of  10  to  20  minutes,  thereby  reducing  the  fluctuation  in  power  by  at  least  20  db  (to  ±3°K).  Fig¬ 
ure  4  provides  an  example  of  a  plot  of  integrated  power  vs  delay.  A  resolution  of  100  psec  (15  km) 
is  achieved  by  making  the  interpulse  period  an  odd  number  of  100-psec  intervals  and  interlacing 
two  time-bases  sampled  at  200-psec  intervals. 

D.  Spectrum  Analyzer 

The  signal  spectra  are  explored  by  a  bank  of  24  filters.  The  outputs  of  the  filters  are 
rectified  with  linear  detectors,  and  these  voltages  are  then  summed  by  carefully  constructed 
analog  integrators.  The  filters  are  driven  by  a  wide-band  (50-kcps)  amplifier  which  is  normally 
switched  off,  and  gated  on  for  a  part  of  the  time-base  equal  to  the  length  of  the  transmitter 
pulse  (0.5  or  1.0msec).  The  delay  t  between  this  gated  portion  and  the  time  at  which  the  trans¬ 
mitter  pulse  is  radiated  determines  the  height  at  which  the  signal  spectrum  is  being  explored. 

In  order  to  calibrate  the  gains  of  the  filter  channels,  the  detectors,  and  to  remove  any  influence 
of  the  earlier  stages  of  the  receiver  on  the  shape  of  the  spectra,,  a  second  gated  portion  of  the 
time-base  is  applied  to  the  filter  bank.  This  second  portion  is  equal  in  width  to  the  first,  but  at 
a  delay  corresponding  to  a  height  from  which  no  detectable  signal  is  expected.  The  outputs  of 
the  detectors  are  switched  to  store  the  voltages  corresponding  to  this  second  gate  in  another  set 
of  24  integrators.  Thus  the  first  set  stores  signal  plus  noise  and  the  second  set  only  noise;  by 
taking  a  ratio  between  the  two  sets  equipment  effects  can  be  eliminated.  Because  the  oscillations 
of  the  filters  must  be  allowed  to  decay  completely  before  the  second  pulse  is  applied,  this  com¬ 
pensation  scheme  imposes  restrictions  on  the  maximum  prf,  which  cannot  be  more  than  50  cps 
for  the  present  filters.  Throughout  1963,  single-pole  filters, were  employed  having  identical 
bandpass  characteristics  with  a  half-power  width  of  500  cps  and  center -to-center  spacing  of 
480  cps. 

Both  halves  of  the  spectrum  should  be  identical.  (Even  if  asymmetry  were  to  exist,  it 
would  be  destroyed  by  the  receiver.)  Thus,  only  one  half  is  normally  measured  and  either  the 
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upper  or  lower  side  band  can  readily  be  selected  for  spectral  analysis.  The  integrators  have  a 
self  time  constant  of  several  hours.  At  the  end  of  a  5 -minute  integration  period,  the  receiver 
pulses  to  the  filter  bank  are  removed  and  the  48  stored  voltages  are  automatically  sampled  in 
turn  and  measured  by  means  of  a  14-bit  digital  voltmeter.  These  voltages,  together  with  the 
time,  are  recorded  on  a  punched  paper  tape. 

IV.  OBSERVING  PROCEDURES 
A.  General 

Observations  were  made  at  intervals  of  about  one  week  throughout  1963.  Each  period  of 

$ 

observation  usually  commenced  near  0900  and  ended  at  about  1700  on  the  following  day.  In 
this  way  the  daytime  behavior  was  observed  twice  during  any  one  run.  Typically,  in  any  one 
hour  the  equipment  was  operated  as  shown  in  Table  II. 


TABLE  II 

TYPICAL  EQUIPMENT  OPERATION 

SEQUENCE 

Pulse  Length 

Computer  Integration 

Spectrum  Delay 

(msec) 

(min. ) 

(msec) 

0.  1 

10 

- 

0.5 

15 

1.5,  2.0,  2.5 

1.0 

15 

3.0,  4.0,  5.0 

0 

5 

(and  6. 0) 

The  sums  of  echo  power  formed  in  the  computer  (e.g.,  Fig.  4)  were  recorded  on  magnetic 
tape  for  later  analysis.  Spectrum  measurements  at  any  one  delay  (Table  II)  normally  occupy 
5  minutes,  thus  to  explore  the  full  range  of  heights  takes  almost  an  hour.  During  the  daytime 
in  summer,  useful  measurements  were  obtained  with  a  delay  of  6  msec,  but  in  winter  at  nighttime 
there  was  insufficient  signal  strength  to  make  measurements  beyond  about  3  msec.  At  the  end 
of  each  hour,  the  RF  drive  to  the  transmitter  was  removed  and  the  measurements  repeated  as 
a  check  on  the  proper  operation  of  the  equipment  in  the  absence  of  a  signal.  The  voltage  outputs 
of  the  integrators  were  recorded  in  binary  fashion  on  a  punched  paper  tape  for  later  analysis. 

It  is  also  possible  to  record  the  receiver  IF  signals  for  spectrum  analysis  at  a  later  time;  this 
can  be  useful  when  it  is  required  to  determine  the  behavior  more  rapidly  than  once  per  hour. 

If  the  electron-to-ion  temperature  T^/T^  is  known  (from  the  signal  spectra),  the  electron 
density  at  that  height  can  be  determined  from  absolute  power  measurements.  These  require 
that  the  effective  antenna  aperture  Aq,  the  transmitter  power  P^,  and  the  receiver  system  tem¬ 
perature  Tg  be  determined. 

Antenna  aperture  is  measured  by  observing  the  radio  source  in  Cygnus  A  which  requires  a 
repositioning  of  the  feed  so  that  the  beam  is  directed  toward  the  south  by  2°  from  the  zenith. 

The  transmitter  power  is  continuously  monitored  by  means  of  power  meters  which  are  connected 
to  couplers  inserted  into  the  waveguide.  The  system  temperature  Tg  is  determined  in  the 
following  way.  A  gas  discharge  tube  providing  a  10,000°K  noise  source  is  coupled  to  the 


*Time  periods  are  EST  unless  otherwise  indicated. 
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receiver  via  a  20 -db  coupler  and  ignited  for  2  msec  at  the  end  of  each  sweep  of  the  time-base. 

A  100 °K  rise  in  system  temperature  can  then  be  observed  in  the  plot  of  integrated  receiver 
output  power  vs  delay  formed  in  the  computer  (Fig.  4).  If  the  20 -db  coupler  is  placed  before  the 
TR  cavity,  any  deterioration  in  the  performance  of  the  TR  tube  can  be  recognized. 

Accuracy  in  determining  the  antenna  efficiency  and  system  temperature  is  governed  chiefly 
by  the  uncertainty  in  the  noise  source  used  as  a  reference.  Several  such  noise  sources  have 
been  employed  but  are  usually  observed  to  differ  by  the  order  of  ±1  db.  A  comparable  error 
may  exist  in  the  determination  of  the  transmitter  power  P  .  Hence,  for  the  most  part,  the 
absolute  value  of  the  echo  power  is  not  used  to  determine  the  electron  density;  instead,  the 
profiles  are  assigned  an  absolute  scale  from  ionosonde  measurements  in  which  fQF2  is  deter¬ 
mined.  The  ionosonde  is  located  within  1km  of  the  backscatter  radar. 

B.  Choice  of  Pulse  Length  (Profile  Measurements) 

The  echo  power  profile  represents  the  convolution  of  the  shape  of  the  pulse  Ah  =  cr/2  with 
the  true  electron-density  profile  N(h).  In  order  to  minimize  the  distortion  of  the  profile  intro¬ 
duced  by  the  pulse,  the  shortest  value  of  pulse  length  r  compatible  with  achieving  sufficient 
echo  power  [Eq.  (2)]  is  used.  For  the  Millstone  Hill  apparatus  this  is  100  ^xsec.  The  thickness 
of  the  F-layer  between  points  of  half-maximum  density  is  of  the  order  of  100  km  or  more  (Fig.  5). 
The  resolution  afforded  by  100-^isec  pulses  (15  km)  is  adequate  to  explore  the  main  features  of 
the  profile,  but  features  having  small  height  extent  (e.g.,  the  FI  ledge)  are  lost. 

During  the  0.5-  and  1-msec  pulse  measurements  the  region  above  400-km  altitude  is  ex¬ 
plored  with  better  signal-to-noise  ratio  (though  poorer  resolution).  Fortunately,  in  these  re¬ 
gions  the  scale  height  is  comparable  to  the  pulse  length,  and  serious  distortion  of  the  profile  is 
not  introduced. 

C.  Choice  of  Pulse  Length  (Spectrum  Measurements) 

The  height  resolution  achieved  in  the  spectrum  measurements  is  considerably  poorer  than 
that  for  the  density  measurements  as  a  consequence  of  the  need  to  employ  long  pulses  in  order 
to  minimize  the  distortion  of  the  spectra.  For  0.5 -msec  pulses  the  height  resolution  is  approx¬ 
imately  75  km,  and  the  width  of  the  pulse  (2kcps)  is  of  the  order  of  one -fifth  the  width  of  the 
spectra  (lOkcps).  If  a  higher  operating  frequency  were  employed,  signal  spectra  would  be 
proportionally  wider  and  shorter  pulses  could  be  used.  However,  in  this  case,  the  ratio 
X/47tX^  would  approach  unity  at  a  lower  height  and  thereby  restrict  the  vertical  extent  of  the 
measurements.  One  solution  of  this  dilemma  seems  to  be  to  employ  an  antenna  directed  at  a 
low  elevation  angle  and  to  restrict  the  vertical  extent  of  the  pulse  by  the  obliquity  of  the  path. 

For  a  pulse  of  length  r  sec  and  a  spectrum  delay  t  sec  between  the  leading  edge  of  the 
transmitter  pulse  and  the  receiver  gate  (of  width  r  sec),  the  lowest  height  from  which  echoes 
are  received  is  (c/2)  (t  —  r)  km,  and  the  uppermost  height  is  (c/2)  (t  +  r)  km.  The  combined 
action  of  the  finite  width  of  the  pulse  and  the  equal  width  of  the  receiver  gate  is  to  ’’weight"  the 
electrons  distributed  between  these  heights  in  a  triangular  fashion,  those  at  (c/2)t  km  affecting 
the  spectrum  the  most.  Actually,  it  is  not  the  electron-density  distribution  but  the  echo  power 
which  is  given  this  weighting  by  the  receiver  gate.  As  is  decreasing  with  delay  t  (for 
most  delays),  the  effective  center  of  the  pulse  is  somewhat  lower  than  ct/2.  Allowance  is  made 
for  this  by  computing  the  effective  height  of  the  pulse  as  a  function  of  delay  t,  from  the  observed 
echo  power  Pf,  vs  height  measurements. 
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V.  DATA  ANALYSIS 


A.  Profile  Measurements 

Plots  of  echo  power  vs  delay  (Fig.  4)  are  recorded  on  magnetic  tape  and  later  analyzed  by 
the  computer  in  the  following  way.  First  the  mean  noise  level  is  established  from  the  values  of 
about  20  points  near  the  center  of  the  time-base,  and  this  mean  is  then  subtracted  from  all  the 
points.  Mean  power  of  the  noise  calibration  pulse  (Fig.  4)  is  next  established  and  used  to  define 
the  absolute  level  of  the  echo  power.  Echo  power  is  then  scaled  up  in  proportion  to  the  square 
of  the  height  to  obtain  an  electron-density  profile  N'(h)  which  is  printed  out  later.  The  computer 
normalizes  the  electron-density  profile  to  the  peak  value  and  plots  it  as  a  function  of  height 
(Fig.  5).  Smooth  curves  are  drawn  through  these  plots  by  hand,  and  a  single  density  profile  is 

extracted  from  the  plots  of  the  0.1-,  0.5-  and  1.0-mSec  pulse  observations  in  the  same  hour. 

46 

Figure  6  presents  an  example  of  such  a  combined  profile  and  has  been  published  previously. 

At  this  point  the  profile  Nf(h)  may  differ  from  the  true  profile  N(h)  if  a  is  a  function  of  height 
but,  as  Fig.  6  shows,  the  error  is  not  extremely  large.  Thus,  by  assigning  the  value  for  fQF2 
observed  on  the  ionosonde  to  the  point  of  peak  electron  density,  an  absolute  scale  is  established 
which  enables  the  approximate  value  of  the  plasma  frequency  f  at  all  other  heights  to  be  esti¬ 
mated.  This  is  necessary  in  order  to  interpret  the  spectrum  measurements  accurately. 

Later,  when  the  values  for  T  /T^  have  been  obtained  as  a  function  of  height  (Fig.  7),  a  first- 
order  correction  can  be  applied  using  Eq.  (4)  to  obtain  the  true  density  profile  N(h).  Figure  8 
shows  the  corrected  version  of  the  profile  shown  in  Fig.  6  after  adjustment,  to  allow  for  the  tem¬ 
perature  inequality  shown  in  Fig.  7. 

B.  Spectrum  Measurements 

The  punched  paper  tape  on  which  the  48  integrator  voltages  have  been  recorded  is  later 
analyzed  by  the  computer.  All  the  values  are  squared  and  a  ratio  is  then  taken  of  each  of  the 
values  of  the  24  signal-plus-noise  to  the  24  noise  powers  to  obtain  the  signal-to-noise  power  at 
each  frequency.  By  subtracting  one  from  these  ratios  the  signal  spectrum  is  obtained.  Cor¬ 
rections  are  applied  to  remove  systematic  differences  between  the  integrators  in  each  pair. 
These  corrections  are  obtained  by  taking  a  mean  of  all  the  ratios  observed  during  the  runs 
when  the  transmitter  is  not  radiating  (Sec.  IV). 

In  order  to  analyze  the  spectrum  measurements  it  is  simply  necessary  to  compare  each 
measured  spectrum  with  a  set  of  theoretical  spectra  to  find  which  agrees  best.  Because  the 
electron-to-ion  temperature  ratio  T  /T.  controls  the  ratio  x  between  the  power  in  the  wing  to 
that  at  the  center  frequency,  this  quantity  (x)  is  scaled  from  the  records.  A  second  quantity  f, 
the  frequency  difference  between  the  center  and  a  point  of  half  peak  power,  is  also  scaled  and 
for  a  given  ion  mass  m.  may  be  regarded  as  a  measure  of  the  ion  temperature  T„  The  effect 
of  the  transmitter  pulse  and  the  frequency  resolution  of  the  receiver  spectrum  analyzer  must 
be  allowed  for,  and  in  order  to  do  this  it  has  been  necessary  to  compute  a  large  number  of 
theoretical  spectra  (e.g..  Fig.  1).  These  spectra  were  then  convolved,  first  with  the  spectral 
distribution  of  the  power  in  the  transmitter  pulse,  and  second  with  the  power-vs -frequency 
response  of  one  of  the  filters  in  the  receiver.  The  effect  of  the  convolution  on  theoretical 
spectra  corresponding  to  given  values  of  T.  and  T^/T^  is  to  lower  x  and  increase  f  as  shown 
in  Fig.  9. 
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A  third  effect,  which  has  not  so  far  been  taken  into  consideration,  is  the  action  of  the  re¬ 
ceiver  gate  on  the  shape  of  the  spectra.  As  mentioned  previously,  the  receiver  gate  is  neces¬ 
sary  to  define  the  height  interval  from  which  echoes  are  received.  The  effective  height  of  the 
pulse  center  is  defined  by  a  weighted  distribution  of  echo  power  within  this  gate,  and  this  quantity 
is  calculated.  However,  the  distortion  introduced  by  the  receiver  gate  in  the  shape  of  the  signal 
spectra  was  not  included  in  the  initial  theoretical  calculations.  From  more  recent  computations, 
it  appears  that  when  the  gate  width  is  made  equal  to  the  pulse  length  the  spectrum  shape  is 
changed  only  by  a  further  lowering  of  the  wing,  and  hence  the  ratio  is  reduced.  It  seems,  there¬ 
fore,  that  the  neglect  of  this  effect  will  not  in  any  way  influence  the  values  of  T^  reported  herein, 
but  Te  and  Te/T^  may  be  underestimated  by  10  to  20  percent.  This  is  probably  comparable  with 
the  experimental  error  at  most  heights  up  to  500  km  and  less  than  the  experimental  error  above 
500km. 


vations  at  a  height  where  the  plasma  frequency  f  is  close  to  5  Mcps 
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The  solid  lines  in  Fig.  9  represent  a  spectrum  analysis  chart  used  for  1-msec  pulse  obser- 

The  need  to  specify  the 

plasma  frequency  stems  from  the  fact  that,  as  the  radio  wavelength  X  approaches  (X^  = 

'kTe/47rNe2  '),  the  spectrum  shape  changes  rapidly.  This  is  illustrated  in  Fig.  10  where  spectrum 
charts  for  a  1-msec  pulse  at  2.5-  and  1.0-Mcps  plasma  frequency  f^  have  been  superimposed. 

Thus,  at  the  wavelength  (68  cm)  employed  here,  the  shape  of  the  spectrum  depends  upon  Tg,  T^ 
and  the  density  N  (which  then  specifies  the  Debye  length).  This  accounts  for  the  need  to  obtain 
an  approximate  profile  before  performing  the  spectrum  analysis  reduction.  Charts  such  as 
Figs.  9  and  10  have  been  prepared  for  1.0-msec  pulses  for  f^  =  1.0,  1.5,  2.0  and  5.0  Mcps,  and 
for  0.5-msec  pulses  at  2.5  and  5.0  Mcps.  Above  f^  =  5.0  Mcps,  X  »  4ttXj^  and  no  further  changes 

in  the  spectrum  shapes  occur.  This  was  tested  by  computing  the  spectrum  shapes  for  f  =  10  Mcps. 

I-  ^ 

All  charts  discussed  so  far  assume  that  O  is  the  only  ion  present. 

The  values  of  T^/T^  corresponding  to  the  measurements  shown  in  Fig.  6  are  shown  in  Fig.  7, 

and  the  separate  dependence  of  T  and  T.  on  this  occasion  is  shown  in  Fig.  11. 

e  1  + 

It  is  possible  that  above  500  km  appreciable  amounts  of  He  were  present  at  the  time  these 
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measurements  were  made.  The  distribution  of  He  ions  has  been  discussed  on  theoretical 
grounds,^'  ^  but  unfortunately  there  is  inadequate  experimental  evidence  at  this  time  to  verify 
the  accuracy  of  these  models.  What  evidence  is  available  is  somewhat  contradictory  as 
Gringauz,  et  al.^  find  a  complete  transition  from  0+  to  He+  over  the  height  range  520  to  620  km. 
The  theory^  and  the  experimental  results  of  Taylor,  et  al.^  indicate  a  more  gradual  transition. 

At  no  time  during  the  course  of  the  backscatter  measurements  have  spectra  been  obtained  that 
definitely  indicate  the  presence,  of  He+  ions,  though  the  spectra  obtained  for  the  uppermost 
heights  (750  km  by  day,  500  km  by  night)  are  uncertain  due  to  the  low  signal-to-noise  ratio,  and 
perhaps  as  much  as  30  percent  He+  might  be  undetected  (Fig.  2).  Figure  12  shows  that  the 
measurements  made  on  2  July  at  720-km  equivalent  height  are  consistent  with  theoretical 
spectra  computed  for  a  plasma  containing  no  He+  (i.e.,  100  percent  0+)  or  one  containing  as 
much  as  20  percent  He+.  Unfortunately,  as  Fig.  13  shows,  the  presence  of  even  small  amounts 
of  He+  requires  a  significant  change  in  the  interpretation  of  the  results.  In  the  specific  example 
shown  in  Fig.  12,  it  is  not  possible  to  distinguish  experimentally  between  two  extreme  cases: 
namely,  (1)  wholly  0+,  T  /T.  =  1.14,  Tg  =  2320°K,  T.  =  2040°K;  and  (2)  80  percent  0+,  20  per- 
cent  He,  T^/T^  =  1.4,  Tg  =  1960°K,  T^  =  1410°K.  These  diverging  interpretations  are  indicated 
in  Fig.  11. 
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On  a  theoretical  analysis  of  the  behavior  of  the  electron  and  ion  temperatures  in  the  upper 
F-region,  it  is  possible  to  decide  between  the  two  possible  curves  shown  in  Fig.  11.  Geisler  and 
Bowhill^0  who  have  performed  such  an  analysis,  show  that  by  far  the  most  important  parameter 
is  the  very  high  thermal  conductivity  of  the  electron  gas  above  about  400  to  500  km.  Thus,  Te 
is  expected  to  be  isothermal  above  this  height  almost  irrespective  of  the  behavior  of  T.  or  the 
neutral  temperature  T^.  If  this  is  correct,  the  solid  curve  of  Fig.  11  (100  percent  0+)  would 
appear  to  be  closest  to  the  true  situation,  and  the  restriction  of  analyzing  the  results  on  the 
assumption  that  0+  predominates  may  be  deemed  to  be  in  order,  where  the  results  show  that 
is  isothermal  at  the  greatest  heights. 

In  principle,  it  should  be  possible  to  determine  the  ratio  0+/Ue+  by  careful  measurement 
of  the  slope  of  the  spectrum  near  the  point  of  half  peak  intensity  (Fig.  2).  In  the  measurements 
described  here,  this  slope  has  been  set  largely  by  the  distortion  introduced  by  the  transmitter 
pulse  and  receiver  filter  characteristics.  It  is  hoped  that  by  using  improved  filters  in  the  re¬ 
ceiver  spectrum  analyzer  this  problem  can  be  minimized. 

In  summary  it  may  be  said  that  the  interpretation  of  the  spectra  corresponding  to  heights  of 

more  than  500  km  is  uncertain  at  the  present  time  owing  to  the  likely  presence  of  He+ions.  Despite 

this.  Fig.  11  shows  clearly  that  T.  increases  markedly  with  height  as  distinct  from  the  neutral 

temperature  T  ,  which  is  widely  thought  to  be  independent  of  height  above  300  km  (Ref.  61).  Either 
n  62  63 

this  view  is  incorrect,  or  T.  >  above  300  km  as  Dalgarno  and  Hanson  have  suggested. 

VI.  ACCURACY  OF  RESULTS 
A.  Profile  Measurements 

Systematic  errors  have  occurred  during  the  course  of  the  profile  measurements  due  to 
poor  recovery  of  the  TR  tube.  Unfortunately,  this  condition  is  usually  not  recognized  during 
an  operation  as  it  is  masked  by  the  ionospheric  echo.  When  this  malfunction  has  been  very 
severe  it  may  be  recognized  in  the  analyzed  profiles  by  the  fact  that  (for  1-msec  pulses)  N(h) 
becomes  negative  near  700  km  Many  of  the  measurements  made  in  1963  suffered  because  of 
this  trouble.  In  October,  provision  was  made  to  monitor  the  recovery  of  the  TR  tube  and  the 
first  stages  of  the  receiver  by  switching  the  receiver  off-tune  sufficiently  (50kcps)  to  place  the 
echo  outside  the  receiver  pass  band.  When  not  too  severe,  the  effect  of  the  TR  tube  trouble  has 
been  to  cause  a  systematic  underestimation  of  N(h)  which  increases  with  height.  An  example 
of  this  behavior  is  given  in  Fig.  14  in  which  a  backscatter  profile  obtained  on  12  July  is  compared 
with  a  density  curve  obtained  from  the  results  of  the  topside  sounder  Alouette  (Van  Zandt,  pri¬ 
vate  communication).  It  can  be  seen  that  serious  error  is  encountered  above  500-km  height. 

Unfortunately,  the  fact  that  this  was  a  continuing  difficulty  was  not  recognized  until  October. 
Earlier  attempts  to  effect  a  cure  were  always  short-lived  as  a  result  of  changes  in  the  charac¬ 
teristics  of  the  discharge  tubes  as  they  aged.  Since  October,  it  is  believed  that  the  density  pro¬ 
files  are  of  comparable  accuracy  to  those  shown  in  Fig.  8.  However,  prior  to  October  it  was  not 
always  possible  to  determine  whether  or  not  the  results  for  a  given  day  were  accurate  or  not, 
unless  the  malfunction  was  particularly  severe.  Thus  the  values  for  the  density  above  500  km 
by  day  and  about  350  km  at  night  should  be  accepted  with  some  caution. 

The  random  errors  of  measurement  increase  rapidly  with  height  (Fig.  5)  because:  (1)  at  the 
uppermost  heights  the  echo  power  is  obtained  as  the  difference  between  two  large  uncertain 
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numbers,  and  (2)  in  computing  the  profile  N’(h)  these  differences  are  scaled  in  proportion  to 
2 

(h)  .  The  profile  measurements  depend  upon  the  determination  of  the  total  echo  power  in  a  wide 
filter,  thus  they  can  always  be  conducted  to  heights  at  which  the  spectrum  measurements  are 
unusable.  Thus,  when  the  equipment  is  operating  properly,  the  upward  extent  of  the  profile 
measurements  is  curtailed  by  the  absence  of  good  spectrum  results  and  not  by  the  errors  in  the 
profiles  themselves.  The  agreement  between  the  rocket  and  backscatter  measurements  shown 
in  Fig.  8  suggests  that  on  this  occasion  there  were  no  serious  sources  of  error  in  the  profile 
measurements.  However,  we  note  that  if  the  uppermost  spectra  (Fig.  12)  are  interpreted  as 
indicating  the  presence  of  20  percent  He+,  the  ratio  T^/T^  would  be  raised  to  1:4  and,  as  a  re¬ 
sult,  the  backscatter  curve  would  be  adjusted  to  agree  completely  with  the  rocket  results  (full 
line)  at  720  km.  Evidently  the  interpretation  of  the  spectra  can  constitute  one  of  the  major 
sources  of  error  in  profile  measurements. 

Below  the  peak  of  the  F-region,  ground  clutter  echoes  restrict  the  lowest  measurements  to 
about  200  km.  These  echoes  are  caused  by  horizontal  radiation  from  the  antenna  feedhorn.  It 
is  possible  that  steps  could  be  taken  to  reduce  these  echoes,  but  the  resolution  obtainable  in  the 
spectrum  measurements  would  be  inadequate  to  explore  this  region,  as  T^/T^,  T^,  and  the  ion 
mixture  all  change  rapidly  with  height.  It  is  planned  to  use  an  oblique  incidence  radar  operating 
at  1295 Mcps  to  explore  the  region  100-  to  300-km  height. 

B.  Spectrum  Measurements 

Each  spectrum  measurement  occupies  a  period  of  five  minutes  during  which  time  15,000 
pulses  are  averaged.  When  the  signals  are  stronger  than  the  noise  (at  low  heights),  the  error 
in  each  point  is  set  simply  by  the  uncertainty  in  the  signal  power  corresponding  to  this  number 
of  samples  (±1  percent).  When  the  signal  is  less  than  the  noise,  the  signal-to-noise  power 
ratio  will  be  times  the  uncertainty  in  the  noise  level  (i.e.,  ±^y~2  percent).  When  the  spectra 
are  later  scaled,  the  determination  of  the  ratio  x  is  always  more  difficult  than  the  half  width 
f,  since  this  is  a  ratio  between  two  like  numbers.  The  error  in  x  might  be  expected  to  be 
about  n/" 2  times  that  in  a  single  point  (i.e.,  ±2  percent  noise/signal  power  ratio  at  low  signal 
levels) . 

In  practice,  when  the  signal  power  is  weak,  systematic  errors  are  introduced  by  low-level 
interference  of  various  forms.  Airborne  altimeters  and  search  radars  appear  to  be  the  worst 
offenders.  In  the  profile  measurements  these  interfering  signals  raise  the  absolute  noise  level 
but  do  not  systematically  distort  the  profiles,  as  they  are  not  synchronized  to  our  repetition 
frequency.  However,  in  the  spectrum  measurements  they  will  reduce  the  apparent  signal-to- 
noise  ratio  over  a  limited  frequency  range.  These  interfering  signals  appear  at  different  fre¬ 
quencies  in  an  unpredictable  fashion  and,  as  a  rule,  useful  spectra  are  not  obtained  when  the 
signal-to-noise  power  ratio  is  less  than  0.2. 

We  have  already  discussed  the  effect  of  an  ion  mixture  on  the  interpretation  of  the  spectra. 
The  poorest  spectra  corresponding  to  the  uppermost  heights  are  in  fact  those  for  which  this 
problem  arises.  At  present,  the  interpretation  of  the  spectra  above  500 -km  height  must  be 
regarded  as  uncertain. 


12 


VII.  OBSERVATIONS  IN  1963 


A.  General 

This  report  is  intended  to  give  a  full  account  of  the  Millstone  Hill  ionospheric  radar  and  the 
backscatter  results  obtained  in  1963.  The  results  presented  in  Sec.  V  have  been  published  pre¬ 
viously,  together  with  other  measurements  made  in  July,  to  demonstrate  the  utility  of  the 
method.4^ 

Table  III  lists  the  dates  and  times  when  observations  were  made.  From  28  February  on¬ 
ward,  the  receiver  was  operated  with  the  synchronously  pumped  electron-beam  parametric 
amplifier  as  first  stage.  No  other  major  modifications  to  the  equipment  were  made  during  the 
report  period,  although  numerous  minor  changes  were  incorporated  to  improve  reliability  and 
performance.  Radar  equipment  failures  occurred  on' few  occasions  and  were  the  least  trouble¬ 
some  cause  of  lost  observing  time. 

Combinations  of  human  and  equipmental  errors  in  operation  of  the  CG  24  computer  caused 
all  profile  results  for  eight  consecutive  weeks  in  May  and  June  to  be  destroyed  and,  in  addition, 
some  shorter  periods  later  in  the  year  were  also  lost  (see  Table  III).  The  next  most  serious 
trouble  was  caused  by  various  kinds  of  external  radio  interference,  which  resulted  in  the  loss 
of  some  observation  time  and  in  the  reduction  of  the  accuracy  of  results  obtained  during  much 
of  the  remaining  time. 

B.  Reduction  Procedure 

The  data  gathered  in  this  program  were  analyzed  according  to  the  methods  outlined  in 
Sec.  V.  A  complete  electron -density  profile  is  obtained  together  with  the  height  dependence  of 
T  and  T^  approximately  every  hour  of  observing  time.  In  order  to  reduce  this  large  amount 
of  data  to  manageable  proportions,  the  electron-density  plots  are  presented  as  diagrams  in 
which  contours  of  constant  critical  frequency  are  drawn  as  functions  of  altitude  and  time. 

One  difficulty  in  presenting  the  diurnal  variation  of  the  temperature  is  that  the  equivalent 
heights  of  the  spectra  vary  throughout  the  day  as  a  consequence  of  changes  in  the  shape  of  the 
layer.  Hence,  any  given  delay  does  not  correspond  to  a  given  fixed  height.  Thus  it  has  proved 
convenient  to  present  the  temperature  results  as  contours  of  constant  T^/T^,  T^  and  T^  on  dia¬ 
grams  having  height  as  ordinate  and  time  as  abscissa. 

In  principle  it  would  be  possible  to  prepare  these  diagrams  for  each  day,  but  in  practice 
the  large  gaps  introduced  by  loss  of  data  through  one  cause  or  another  make  it  preferable  to 
average  many  days'  data  into  a  single  plot.  The  diurnal  and  seasonal  behavior  were  considered 
the  most  important  characteristics  to  study.  Consequently,  the  results  have  been  used  to  con¬ 
struct  mean  plots  showing  the  behavior  averaged  over  each  calendar  month.  In  order  to  do  this, 
the  day  is  treated  as  24  separate  hourly  intervals,  and  a  mean  taken  of  all  the  results  obtained 
in  each  given  hour  over  a  given  month.  In  the  case  of  the  electron-density  profiles  it  is  first 
necessary  to  adjust  all  the  profiles  to  their  mean  height  before  computing  the  mean  shape.  The 
mean  profile  is  then  accorded  the  average  value  of  the  density  at  its  peak.  When  the  day-to-day 
critical  frequency  is  subject  to  rapid  change  (e.g.,  during  the  equinoctial  periods),  some  difficulty 
can  be  encountered  in  constructing  the  diurnal  plots  if  records  are  missing  for  parts  of  certain 
days.  Thus,  gaps  are  left  on  several  of  these  plots  where  there  are  insufficient  data  to  make 
meaningful  averages. 
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TABLE  111 

IONOSPHERIC  BACKSCATTER  OBSERVATIONS  -  1963,  1964 

MILLSTONE  RADAR 

Sfort  Time 

End  Time 

Start  Time 

End  Time 

Date 

(EST) 

(EST) 

Comments 

Date 

(EST) 

(EST) 

Comments 

1963 

1963 

A/5  Feb 

1100 

1700 

2100-0400  receiver  malfunction,  no  data; 

16/17  Aug 

0630 

1550 

Severe  interference; 

diode  parametric  amplifier  in  use 

NO  PROFILE  DATA  2000-0500 

19/20  Feb 

1100 

1615 

Electron-beam  parametric  amplifier  used 

23/24  Aug 

0900 

1640 

NO  PROFILE  DATA  0100-0600 

28  Feb/1  Mor 

1030 

1600 

Electron-beam  parametric  amplifier 

29/30  Aug 

1100 

1100 

synchronously  pumped 

7/8  Mor 

1230 

1700 

6/7  Sept 

0820 

1540 

14/15  Mor 

1040 

1600 

1300-1700  on  14  Mor,  NO  DATA 

13/14  Sept 

0915 

1600 

Strong  auroral  interference 

(transmitter  malfunction) 

NO  PROFILE  DATA  0000-1600 

20/21  Mar 

0930 

1400 

20  Sept 

0740 

1600 

28  Mar 

0500 

1915 

21/22  Sept 

0640 

0700 

Interference;  NO  PROFILE  DATA  0000-0600 

4/5  Apr 

0930 

1415 

27  Sept 

0740 

1600 

10/11  Apr 

0920 

1550 

28/29  Sept 

0740 

0800 

TR  trouble  severe; 

NO  PROFILE  DATA  2000-0600 

17/18  Apr 

1210 

1815 

4/5  Oct 

1100 

1540 

Digital  voltmeter  foiled;  all  spectrum  dota 

25/26  Apr 

1140 

1315 

NO  PROFILE  RESULTS 

recorded  and  processed  on  29/30  Oct 

7/3  May 

0740 

1200 

NO  PROFILE  RESULTS 

18/19  Oct 

0840 

1530 

NO  PROFILE  DATA  1000-1530  on  19  Oct 

10/11  Moy 

1550 

1800 

NO  PROFILE  RESULTS 

25/26  Oct 

0650 

1530 

Several  equipment  malfunctions; 

NO  PROFILE  DATA  0400-0930 

31  May/1  June 

1015 

1600 

NO  PROFILE  RESULTS 

1/2  Nov 

1130 

1650 

NO  PROFILE  DATA  0100-0530 

7/8  June 

1000 

1540 

NO  PROFILE  RESULTS 

7/8  Nov 

1000 

1700 

NO  PROFILE  DATA  1900-0600 

14/15  June 

0840 

1600 

NO  PROFILE  RESULTS 

14/15  Nov 

1045 

1620 

21/22  June 

0630 

1600 

NO  PROFILE  RESULTS 

21/22  Nav 

0930 

1640 

28/29  June 

0930 

1530 

NO  PROFILE  RESULTS 

until  0600,  29  June 

29/30  Nov 

0910 

1715 

7/8  July 

0900 

1230 

6  Dec 

0930 

1300 

Punch  failure  coused  lass  of  spectrum 

data  obtained  ofter  1 300 

11/12  July 

0600 

1530 

13/14  Dec 

0930 

1630 

NO  PROFILE  DATA  1935-0010 

19  July 

0630 

1940 

20/21  Dec 

1030 

1715 

NO  PROFILE  DATA  0240-0720 

21  July 

0820 

1940 

Severe  interference  0800-1400; 

NO  PROFILE  DATA  0900-1300 

27/28  Dec 

0920 

1645 

26/27  July 

0600 

1600 

1964 

3/4  Jon 

0910 

1600 

NO  SPECTRUM  DATA  0030-1015 

7/3  Aug 

0640 

1600 

TR  trouble  severe; 

NO  PROFILE  DATA  0100-0500 

10/11  Jon 

0930 

1630 

NO  SPECTRUM  DATA  1400-1630  on  11  Jon 

9/10  Aug 

0610 

1620 

Severe  interference; 

17/18  Jan 

1000 

1650 

NO  SPECTRUM  DATA  1000-1500  an  17  Jon 

NO  PROFILE  DATA  2000-2300 
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No  large  changes  in  the  height  behavior  of  the  ratio  Te/T.  were  encountered  as  a  function 
of  date;  hence,  some  saving  in  reduction  time  was  achieved  by  taking  a  mean  of  the  uncorrected 
N'(h)  profiles  and  correcting  this  using  the  mean  Te/T.  curve  obtained  for  that  month. 

The  quality  of  the  spectrum  results  varied  markedly  with  altitude,  time  of  day,  and  depend¬ 
ing  on  whether  or  not  any  radio  interference  was  present.  This  is  accounted  for  in  the  compu¬ 
tations  of  the  mean  temperatures  by  weighting  the  values  of  Tg,  T^  and  T  /T.,  one,  two,  or 
three  times  according  to  the  confidence  that  could  be  placed  in  the  corresponding  spectrum. 

This  weighting  is  somewhat  arbitrary,  but  does  serve  to  bring  the  mean  values  for  T^  and  T. 
closer  to  the  values  obtained  from  the  best  spectra.  Thus  a  weighted  mean  is  obtained  for  the 
values  of  T^,  T^/T^,  and  T.  obtained  in  each  hour,  in  each  month,  for  each  given  delay.  This 
mean  is  assigned  the  weighted  mean  equivalent  height.  For  each  hour  a  plot  is  then  constructed 
of  Tg  and  T^  as  a  function  of  height,  and  smooth  curves  are  drawn  through  the  points  because 
no  sharp  discontinuities  in  either  of  these  temperatures  would  be  physically  plausible.  In  this 
way  the  effects  of  single  spurious  points  are  minimized  unless  they  happen  to  be  uppermost.  In 
constructing  these  curves,  greatest  weight  is  given  to  the  points  representing  the  average  of  the 
largest  number  of  determinations.  Finally,  a  new  determination  of  the  ratio  T^/T.  with  altitude 
is  obtained  from  these  smooth  T^  and  T.  curves  and  compared  with  the  weighted  mean  values. 
Usually  good  agreement  would  be  found,  but  where  there  is  a  difference  the  Te/T^  curve  derived 
from  the  two  temperature  curves  is  followed  instead  of  the  original  points. 

No  effort  was  made  during  the  data-taking  operations  to  synchronize  the  measurements  to 
fixed  times;  hence,  it  is  assumed  that  the  average  values  are  representative  of  conditions  at  the 
midpoint  of  each  hour,  and  they  are  plotted  accordingly. 

Where  the  electron-density  profiles  were  destroyed,  or  otherwise  missing,  the  spectrum 
measurements  could  still  be  analyzed  if  the  approximate  critical  frequency  at  the  corresponding 
height  could  be  estimated.  To  do  this,  it  was  assumed  that  the  plasma  frequency  was  the  same, 
on  the  day  for  which  it  was  missing,  as  actually  observed  on  the  nearest  date  at  that  hour.  It 
is  believed  that  the  errors  inherent  in  this  procedure  are  not  sufficient  to  cause  serious  inac¬ 
curacy  in  interpretation  of  the  spectra. 

VIII.  ELECTRON-DENSITY  MEASUREMENTS 

Figures  15(a)  through  (j)  show  results  of  the  electron-density  determinations.  The  months 
of  May  and  June  are  missing  for  the  reasons  mentioned  previously  (see  Sec.  VII).  Other  months 
have  gaps  where  the  paucity  of  data  makes  it  difficult  to  construct  meaningful  average  profiles. 
No  attempt  has  been  made  to  smooth  the  plots  shown  in  Figs.  15(a)  through  (j),  instead  the  values 
for  each  hourly  interval  have  simply  been  joined  by  straight  lines.  The  main  features  of  these 
plots  are  discussed  briefly  here. 

There  is  little  summer-winter  difference  in  the  density  at  the  maximum  of  the  layer  because 

sunspot  minimum  conditions  were  being  approached  at  this  time.  The  highest  midday  critical 

frequency  observed  was  7.0Mcps  in  November,  and  the  lowest  was  5.0Mcps  in  July  and  August. 

Millstone  Hill  is  located  under  the  "winter  spur"  of  high  latitude  ionization  but  this  is  not  pro- 

64 

pounced  near  sunspot  minimum.  The  critical  frequency  rises  rapidly  at  dawn  and  usually 

reaches  a  maximum  near  noon  or  a  little  before.  In  all  months,  hmax  is  lowest  at  about  2  to  3 

hours  after  sunrise  (around  0700  to  0900)  and  then  has  a  value  lying  between  220  and  230  km. 

As  the  day  advances,  h  rises  almost  linearly  to  about  310  km  which  it  reaches  a  little  after 

max 

midnight.  At  sunrise,  h  decreases  rapidly. 

°  max 
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In  summer  and  the  equinoctial  months,  the  shape  of  the  layer  is  roughly  constant  throughout 
the  major  part  of  the  day.  The  rise  in  causes  the  layer  as  a  whole  to  increase  in  height 

until  sunset.  During  the  evening  the  shape  of  the  layer  changes,  but  remains  roughly  constant 
during  the  early  morning  hours  —  when  the  whole  layer  can  be  seen  to  be  disappearing.  In  winter 
(e.g.,  November)  there  is  no  time  when  the  layer  shape  is  not  changing. 

A  prominent  feature  of  the  variation  of  fQF2  on  the  east  coast  of  the  United  States  is  a  large 
increase  occurring  at  sunset.  Often  fQF2  will  then  exceed  its  value  at  any  other  time  of  the  day. 
In  summer  months  (e.g.,  July)  the  noon  and  evening  maxima  may  be  separated  by  as  much  as 
10  hours,  but  in  the  equinoxes  (March  and  October)  this  is  reduced  to  about  6  hours.  The  eve¬ 
ning  increase  is  absent  in  winter.  Cause  of  the  evening  increase  is  very  evident  from  inspection 
of  the  results  for  July  and  August.  At  650-km  altitude,  density  begins  to  decrease  2  to  3  hours 
before  the  evening  increase  reaches  its  peak.  A  similar  decrease  in  density  takes  place  at  lower 
heights,  somewhat  delayed  after  that  at  650  km.  Evidently  ionization  moves  downward  and  piles 
up  at  h  .  The  reason  for  the  downward  motion  of  ionization  at  this  time  is  thought  to.be  the 
large  fall  in  T^  at  sunset,  which  causes  the  equilibrium  scale  height  H  to  decrease  rapidly. 

Similar  behavior  was  observed  during  the  eclipse  of  20  July  1963  (Ref.  47).  Increases  of  f  F2 
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during  eclipses  have  been  shown  to  depend  upon  the  eclipse  being  total  at  FI  region  heights, 

and  the  magnetic  dip  being  large  (I  $.60°).  We  suggest  that  this  latter  condition  is  most  impor¬ 
tant  for  observing  an  evening  increase.  When  I  is  large,  fast  photoelectrons  can  travel  upward 
to  great  heights  during  the  daytime,  and  thus  give  rise  to  high  values  for  the  exospheric  electron 

temperature.  Also,  the  downward  diffusion  (which  depends  upon  the  dip  angle  I  according  to 
2 

sin  I)  can  be  rapid  at  sunset. 

In  a  separate  paper,48  the  appearance  of  the  midlatitude  evening  increase  in  fQF2  will  be 
discussed  in  some  detail.  The  absence  of  the  effect  in  winter  can  be  attributed  to  the  more 
gradual  and  less  extensive  fall  in  Tg  at  sunset.  Also,  less  ionization  is  available  above  hmax 
to  participate  in  downward  diffusion.  At  sunspot  maximum,  the  effect  is  less  pronounced  but 
does  serve  to  maintain  high  values  of  fQF2  almost  until  midnight.  This  change  between  sunspot 
minimum  and  maximum  can  be  attributed  to  reduction  in  the  temperature  difference  T^  —  T^, 

thought  to  occur  at  sunspot  maximum  as  a  consequence  of  the  higher  electron  density 66  The 
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effect  of  the  magnetic  declination  on  the  effect  reported  by  Eyfrig,  *  if  real,  must  now  be 
interpreted  as  evidence  for  electrodynamic  drifts  which  either  oppose  or  assist  the  downward 
diffusion  induced  by  the  temperature  change. 

Systematic  increases  in  fQF2  are  observed  at  night  during  winter  months  (November,  De¬ 
cember  and  January).  These  increases  all  begin  between  midnight  and  0200,  well  before  local 
sunrise.  There  is  some  evidence  for  an  increase  in  fQF2  in  March  beginning  at  0400  —  again 
before  sunrise.  We  can  exclude  the  possibility  that  this  is  a  consequence  of  freshly  created 
ionization,  perhaps  produced  by  cosmic  rays  precipitating  into  the  atmosphere.  If  this  were 
the  case  we  should  expect  a  pronounced  increase  in  T^  at  this  time.  As  we  shall  see,  T^  usually 
decreases  during  these  periods,  and  this  can  be  explained  as  a  result  of  the  inability  of  the  night¬ 
time  heat  source  (to  be  discussed)  to  maintain  the  same  temperature  in  the  presence  of  additional 
ionization. 

It  also  seems  unlikely  that  this  phenomenon  is  identical  in  nature  to  the  summer  evening  in¬ 
crease.  On  the  basis  of  such  a  model,  it  would  be  exceedingly  difficult  to  explain  why  the 
increase  occurs  in  the  early  hours  of  the  morning,  and  why  large  decreases  in  T^  and  T^  which 
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could  give  rise  to  such  an  effect  are  not  seen.  Nor  can  the  increase  in  fQF2  be  attributed  simply 
to  a  large  change  in  the  shape  of  the  layer,  which  causes  a  redistribution  of  ionization  in  which 
more  electrons  appear  at  the  peak.  The  increase  observed  at  the  peak  is  evident  at  other  levels 
above  and  below  h 

max 

The  regular  evening  decrease  in  density  appears  to  be  arrested  near  2300  in  these  winter 

months,  and  the  density  at  all  levels  then  increases  reaching  peaks  at  about  0230  in  the  400-  to 

500 -km  region  and  at  about  0400  at  h max ♦  This  suggests  an  influx  of  ionization  from  great 

heights,  which  must  commence  near  2300  and  persist  for  several  hours.  The  phenomenon  is 

common  at  midlatitude  stations  in  winter  at  sunspot  minimum.  It  so  happens  that  the  evening 

increase  discussed  above  is  most  pronounced  at  these  same  stations  in  summer  at  sunspot 

minimum.  We  suggest  that  the  early  morning  winter  increase  is  a  consequence  of  the  late 

evening  increase  that  occurs  in  the  summer  (southern)  hemisphere.  Without  needing  to  invoke 

69 

actual  transport  of  ionization  from  one  hemisphere  to  another  as  proposed  by  Rothwell,  a 
plausible  explanation  (though  highly  speculative)  can  be  offered.  We  suppose  that  above  about 
400  km  the  electron  temperature  along  a  field  line  is  approximately  constant  as  a  result  of  the 
high  thermal  conductivity  of  the  electrons.  During  the  northern  winter  evening  hours,  a  large 
region  of  the  exosphere  must  be  warmed  from  the  southern  foot  of  the  field  line  where  the  sun 
has  not  set.  When  the  sun  does  set,  the  temperature  along  the  whole  field  line  may  decrease 
giving  rise  almost  immediately  to  a  pronounced  increase  in  fQF2  in  the  summer  hemisphere 
and  to  a  considerably  smaller  increase  in  the  winter  hemisphere  at  some  later  time.  This  idea 
has  been  discussed  in  more  detail  elsewhere.70 

The  total  electron  content  of  the  ionosphere  can  be  determined  by  Faraday  rotation  measure - 
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ments.  The  two-frequency  moon-reflection  measurements  are  probably  most  accurate.  At 

temperate  latitudes  in  winter,  the  daytime  ratio  na/n^  of  the  number  n&  of  electrons  above 

h  F2  to  the  number  below  n,  is  approximately  2.5  to  1,  and  at  sunspot  minimum  the  summer 

day  ratio  is  about  twice  this.  We  are  not  in  a  position  to  test  these  values  as  the  profiles  do 

not  extend  to  either  sufficiently  low  or  sufficiently  high  altitudes.  However,  it  is  clear  that, 

although  fQF2  decreases  from  winter  to  summer,  the  density  at  a  height  of  say  600  km  increases 

from  winter  to  summer.  Thus  the  scale  height  H  of  the  ionization  above  hmaxF2  is  larger  in 

summer  than  winter,  and  this  causes  the  increase  in  the  ratio  n  /n,  . 

a  b 

IX.  ION-TEMPERATURE  MEASUREMENTS 

Plots  were  constructed  by  drawing  smooth  temperature  curves  for  each  hourly  interval  and 
plotting  as  points  the  height  at  which  a  given  temperature  occurred;  these  points  were  then 
joined  with  straight  lines.  Fluctuations  occur  in  these  plots  as  a  result  of  experimental  errors, 
and  a  better  presentation  can  be  obtained  by  drawing  smooth  curves  to  follow  the  general  be¬ 
havior  of  the  lines  in  these  plots.  These  smoothed  plots  are  shown  in  Figs.  16(a)  through  (1). 

It  was  frequently  difficult  to  decide  to  what  extent  the  original  curves  should  be  smoothed,  and 
it  is  likely  that  some  real  variations  have  been  lost  in  this  process.  Nevertheless,  the  smoothed 
plots  do  retain  the  main  features  of  the  diurnal  behavior.  The  values  for  heights  above  500  km 
are  possibly  too  large  as  a  result  of  the  influence  on  the  spectra  of  He  ions,  and  as  such  must 
be  accepted  with  caution. 

The  principal  features  of  these  plots  are  the  rapid  increase  in  T^  at  dawn,  and  somewhat 
less  rapid  decrease  at  sunset.  These  changes  occurred  at  all  heights  above  about  250  km. 
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During  the  daytime  the  ion  temperatures  show  an  increase  of  about  2°/km  in  the  height  range 
250  to  500  km,  irrespective  of  season.  Above  500  km,  T.  increases  by  more  than  2°/km  in  the 
daytime,  suggesting  that  the  ion  temperature  is  converging  toward  the  electron  temperature. 

At  night,  T\  increases  by  about  l°/km  up  to  about  500  km,  and  usually  by  less  above  this  height. 
There  is  evidence  in  some  of  the  plots  for  sudden  decreases  of  during  the  hours  of  darkness 
(during  the  winter  months),  and  these  appear  to  be  associated  with  increases  in  density  at  that 
time.  Seasonal  variations  of  T.  will  be  discussed  later  (Sec.  XII). 

X.  ELECTRON-TEMPERATURE  MEASUREMENTS 

Figures  17(a)  through  (1)  show  determinations  of  the  electron  temperature.  These  plots 
were  constructed  in  the  same  manner  as  those  for  the  ion  temperature.  It  is  at  once  evident 
that  T^  >  T^  at  all  heights  most  of  the  time.  Like  the  ion  temperature,  T^  shows  a  rapid  in¬ 
crease  at  dawn  and  less  rapid  decrease  at  sunset.  There  is  some  tendency  for  the  electron 
temperature  (at  say  300-km  altitude)  to  reach  its  highest  value  at  about  2  to  3  hours  after  sun¬ 
rise,  but  there  are  no  large  changes  in  T^  at  this  altitude  during  the  major  part  of  the  day.  Un¬ 
like  T.,  the  electron  temperature  increases  rapidly  with  height  (~2°/km)  up  to  about  400  km 

1  + 
and  then  becomes  almost  isothermal  in  the  daytime.  The  influence  of  neglected  He  ions  will 

be  to  overestimate  T^  but  underestimate  Te/T^,  thus  T^  is  likely  to  be  the  quantity  least  in 

error.  Despite  this,  we  feel  the  fact  that  the  region  above  500  km  appears  isothermal  suggests 

that  the  spectra  have  been  properly  interpreted  for  the  most  part;  i.e.,  little  He+  is  present 

during  the  daytime  at  these  altitudes. 

Nighttime  increases  in  appear  to  have  taken  place  during  some  of  the  months.  A  clear 
example  is  that  at  0200  to  0300  in  March.  In  winter  months  the  spectra  were  usually  poor  be¬ 
cause  of  the  low  electron  density.  Partly  as  a  result  of  this,  and  also  because  of  real  fluctua¬ 
tions  in  T^  with  time,  it  has  been  difficult  to  construct  the  nighttime  behavior  in  the  winter 
months  (e.g.,  February).  Precisely  because  the  electron  density  is  low  during  these  months, 
the  electron  temperature  can  be  raised  considerably  by  only  a  small  heat  flux. 

XI.  MEASUREMENTS  OF  ELECTRON-TO-ION  TEMPERATURE  RATIO 

Plots  of  the  ratio  of  T  /T.  are  shown  in  Figs.  18(a)  through  (1).  The  highest  value  of  the 
ratio  is  usually  encountered  at  about  350-  to  400 -km  altitude  shortly  before  noon.  Often  there 
is  a  second  maximum  about  noon  or  later.  In  most  months  these  maxima  have  a  value  of  2.2  or 
2.4,  and  the  highest  ratio  at  this  altitude  was  3.0  in  March.  During  the  daytime,  T^/T^  ^.1.6 
at  all  altitudes  under  observation;  hence,  there  are  no  violent  changes  in  the  ratio  with  altitude. 
The  ratio  T^/T.,  increases  initially  with  altitude  because  Te  is  rising  faster  than  T^,  but  sub¬ 
sequently  decreases  as  Te  becomes  nearly  isothermal  and  T^  begins  to  increase  rapidly.  It  is 
believed  that  these  values  of  T  /T.  are  correct  to  ±20  percent  up  to  about  500  km.  Above  this 
altitude,  the  presence  of  He+  ions  (particularly  at  night)  may  cause  the  spectra  to  be  interpreted 
as  indicating  too  low  a  value  of  Te/T.. 

Although  there  is  a  clear  increase  in  T^/T^  at  sunrise  and  a  corresponding  decrease  at 
sunset,  very  rarely  is  thermal  equilibrium  completely  established  at  night.  Because  T^  and  T^ 
tend  to  vary  during  the  night,  their  ratio  T^/T^  often  provides  a  clearer  indication  of  nighttime 
heating.  Instances  of  such  heating  are  evident  in  March  (near  0300),  in  May  (near  0300),  in 
June  (near  0200),  in  October  (near  midnight),  and  in  January  1964  (near  2100).  Even  at  times 
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other  than  near  these  definite  nighttime  increases,  only  during  the  summer  months  is  Te/T. 
ever  as  low  as  1.2.  A  more  typical  ratio  would  be  1.6,  but  the  marked  fluctuations  of  T^/T^  at 
night,  particularly  during  the  winter  months,  make  it  difficult  to  construct  meaningful  plots. 

The  possibility  of  very  large  values  of  T  /T.  occurring  at  sunrise  at  certain  heights  has 
6>  2  c  i 

been  raised  by  Dalgarno.  Evidence  of  this  during  May  must  be  discounted,  as  there  is  reason 
to  believe  there  was  appreciable  heating  during  the  early  morning  hours  in  this  month.  If  an 
effect  associated  with  sunrise  does  occur,  it  might  be  expected  chiefly  in  the  winter  months 
when  the  predawn  electron  density  is  lowest,  yet  it  does  not  seem  to  have  been  observed.  This 
implies  that  the  local  heating  assumption  invoked  by  Dalgarno  does  not  hold  in  the  ionosphere. 

XII.  AVERAGE  DAYTIME -TEMPERATURE  BEHAVIOR 

In  most  months,  fairly  stable  temperature  conditions  seem  to  exist  throughout  the  major 
part  of  the  day.  Accordingly,  a  period  of  6  or  7  hours  was  selected  during  the  middle  of  the 
day  in  each  month  when  the  temperatures  appeared  to  change  least,  and  mean  temperature  pro¬ 
files  were  constructed  for  these  periods.  These  profiles  are  shown  in  Figs.  19(a)  through  (1) 
where,  in  general,  the  behavior  of  Te  and  T^  with  height  does  not  have  any  very  marked  seasonal 
dependence,  though  some  changes  can  be  observed. 

Seasonal  variations  in  are  shown  in  Fig.  20  where  the  monthly  mean  values  for  300,  400 
and  500  km  extracted  from  Figs.  19(a)  through  (1)  are  plotted.  It  can  be  seen  that  T.  increases 
from  about  960°K  in  midwinter  to  about  1050°K  in  midsummer  at  300  km.  At  400km  the  varia¬ 
tion  is  about  the  same,  i.e.,  from  1020°  to  1160°K,  but  the  winter-to-summer  variation  at  500km 
is  much  larger  —  approximately  240 °K.  At  300 -km  altitude,  the  values  of  T^  obtained  here  are 
probably  close  to  the  neutral  temperature  T^.  However,  satellite  results  indicate  that  T^  is  a 
maximum  in  the  equinoxes73  and  not  in  the  summer  as  found  for  T . .  It  is  possible  that  some  of 
the  fluctuation  of  the  points  in  Fig.  20  is  significant,  but  there  seems  no  way  of  testing  this  at 
the  present  time.  The  possibility  that  the  summer  increase  in  T^  is  a  consequence  of  the  sum¬ 
mer  measurements  being  made  on  magnetically  disturbed  days  has  been  investigated.  The  mag¬ 
netic  character  figure  A^  is  plotted  in  Fig.  21  for  all  the  days  on  which  measurements  were 
made;  also  shown  is  the  monthly  mean.  It  can  be  seen  that  there  is  no  correlation  between  the 
mean  monthly  temperatures  and  the  mean  values  for  A^.  This  probably  signifies  that  the  largest 
part  of  the  variation  in  Fig.  20  is  simply  a  seasonal  one,  and  a  more  careful  test  needs  to  be 
made  to  determine  whether  or  not  there  are  any  changes  during  disturbed  conditions.  Such  a 
test  is  described  in  Sec.  XIV.  Also,  it  can  be  seen  from  Fig.  20  that  the  solar  ultraviolet  flux 
as  characterized  by  the  sun's  radio  radiation  at  10.7  cm  (S1Q  ?)  was  roughly  constant  throughout 
1963,  so  that  the  summer  increase  in  T^  cannot  be  ascribed  to  changes  in  the  intensity  of  the 
solar  ultraviolet.  It  would  seem  that  this  is  simply  a  consequence  of  the  reduced  solar  zenith 
distance  in  summer. 

The  seasonal  dependence  in  the  electron  temperature  is  shown  in  Fig.  22  for  300-  and  500-km 
altitude.  The  scatter  of  the  points  for  the  electron  temperature  is  larger  than  that  for  T^,  mak¬ 
ing  interpretation  more  difficult.  It  is  clear,  however,  that  T^  is  at  a  minimum  in  the  winter 
months  and  at  a  maximum  in  the  late  spring  to  early  summer.  It  is  probable  that  there  is  a 
second  smaller  maximum  in  the  late  summer.  The  difference  between  maximum  and  minimum 
is  ~500°  to  600°.  The  variation  of  the  maximum  value  of  Te/T^  is  also  shown  in  Fig.  22,  and 
it  seems  clear  that  it  is  this  variation  which  gives  rise  to  the  two  maxima  in  T^. 


19 


The  seasonal  variation  in  can  readily  be  accounted  for  by  the  variation  of  the  sun’s  zenith 
distance  provided  it  is  assumed  that  above  300  km.  A  simple  explanation  can  also  be  of¬ 

fered  for  the  variation  in  Te/T^.  At  first  sight  it  would  seem  that  T^/T^  should  peak  in  the  sum¬ 
mer,  because  the  electron  density  is  then  least  at  the  peak  of  the  F2  layer  and  the  fast  photo¬ 
electrons  should  be  able  to  traverse  the  peak  with  less  energy  loss  than  in  the  winter.  However, 
an  opposing  effect  exists,  namely  that  the  thickness  of  the  layer  is  at  maximum  in  summer,74  and 
perhaps  this  more  than  offsets  the  change  in  density  at  the  peak. 


XIII.  AVERAGE  NIGHTTIME-TEMPERATURE  BEHAVIOR 

Daytime  behavior  is  usually  fairly  stable  and,  by  taking  an  average  over  several  hours,  we 
have  been  able  to  identify  certain  seasonal  trends.  Nighttime  behavior,  on  the  other  hand, 
seems  erratic  (particularly  in  the  winter  months);  hence,  averages  have  been  taken  not  simply 
as  a  means  of  reducing  experimental  errors  but  in  order  to  gain  a  better  understanding  of  the 
average  behavior.  These  nighttime  averages  are  shown  in  Figs.  23(a)  through  (1).  Considerable 
variability  exists  from  month  to  month  reflecting  the  random  changes  within  the  month,  and  also 
the  reduced  experimental  accuracy  at  night  when  the  signals  are  weak.  Despite  this,  some 
regularities  can  be  discerned. 

Below  300  km,  T.  and  T^/T^  are  usually  found  to  be  increasing,  but  above  this  height 
usually  increases  by  less  than  l°/km.  Te/T^  does  not  change  by  very  large  amounts  with 
height,  but  may  exhibit  a  maximum  of  350-  to  400-km  altitude.  Because  Tg/T^  does  not  change 
markedly  with  height,  the  curve  for  tends  to  follow  the  curve  for  T„ 

Figure  24  shows  the  variation  of  T^,  T^,  and  T^/T^  at  350-km  altitude.  If  any  nighttime 

heating  agency  is  active  throughout  the  year  without  any  significant  change  in  flux,  it  would  be 

expected  that  T  /T.  would  be  a  minimum  in  summer  when  the  electron  density  N  is  highest 
^  i  2  63 

(since  the  heat  loss  from  electrons  to  ions  proceeds  at  a  rate  proportional  to  N  ).  This  appears 

to  be  the  case  —  the  lowest  values  of  T^/T^  being  observed  in  July.  The  ion  temperature  shows 

a  pronounced  peak  in  September,  and  one  is  reminded  that  this  was  the  month  that  was  most 

magnetically  disturbed  (Fig.  21).  The  electron  temperature  is  also  large  in  this  month,  and  it 

is  tempting  to  conclude  that  there  is  a  heat  source  whose  magnitude  depends  upon  the  magnetic 

character  figure  A^  but  that  the  daytime  temperatures  are  insensitive  to  this  source,  as  solar 

ultraviolet  is  much  larger. 


XIV.  TEMPERATURE  EFFECTS  ASSOCIATED  WITH  MAGNETICALLY 
DISTURBED  CONDITIONS 

We  have  been  concerned  in  this  study  principally  with  determining  the  normal  daily  and  sea¬ 
sonal  behavior  of  the  quantities  measured;  hence,  magnetic  disturbance  effects  have  not  been  ex¬ 
amined  in  detail.  It  is,  however,  important  to  assess  to  what  extent  magnetic  disturbances 
may  have  played  a  part  in  controlling  the  results.  We  have  already  concluded  that  during  the 
daytime  Tg  and  T^  are  insensitive  to  disturbances  as  characterized  by  the  figure  A^.  Also  we 
have  seen  that  Te  and  T^  apparently  increase  on  disturbed  nights,  but  this  conclusion  rests  on 
the  high  values  observed  at  night  in  September  —  the  most-disturbed  month.  To  test  these  con¬ 
clusions  we  have  divided  September  into  the  quietest  days  (A  ^  14)  and  most-disturbed  days 
(Ap  ^  19).  The  average  values  of  A^  for  the  two  groups  were  11  and  67,  respectively.  The 
mean-temperature-vs-height  dependence  for  the  two  groups  is  shown  in  Figs.  25  and  26.  As 
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far  as  can  be  determined,  is  the  same  on  both  the  quiet  and  disturbed  days  in  the  daytime 

though,  below  400  km,  T.  was  about  100°  higher  on  the  disturbed  day.  At  night,  however, 

extremely  large  differences  appear.  Above  350  km,  T  is  raised  over  800°K  and  T.  by  about 

75  e  1 

200°K.  We  have  discussed  elsewhere  the  behavior  of  the  temperature  on  individual  days  in 
this  month,  and  will  here  only  summarize  the  results.  For  the  daytime  period  (1000  to  1500) 
there  seemed  little  significant  change  in  at  any  height,  but  a  decrease  in  on  disturbed 
days  (Table  IV).  At  night,  and  increased  at  all  heights  during  disturbed  conditions.  The 
average  increase  in  T.  per  unit  increase  of  the  index,  AT^/AK^  ~  190°  (Table  IV)  and 
AT^/AKp,  was  twice  as  large.  Also  given  in  Table  IV  are  the  temperature  increases  as  func¬ 
tions  of  the  three-hour  index  a^.  It  seemed  from  the  scatter  of  the  results  that  the  temperature 
was  more  closely  related  to  A^  than  K^,  though  there  were  insufficient  data  to  be  certain  of  this. 

The  values  of  AT  /AK  given  in  Table  IV  presumably  depend  to  a  very  large  extent  upon  the 
e  p 

electron  density,  and  hence  would  be  expected  to  exhibit  both  diurnal  and  seasonal  changes  —  even 

if  Kp  is  a  proper  index  to  employ  as  a  measure  of  the  heat  flux.  For  example,  we  should  expect 

quite  small  heat  fluxes  to  raise  T  considerably  at  night  in  winter  but  not  in  summer. 
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Evans  and  Taylor  observed  from  Faraday  rotation  measurements  that  during  the  night  in 

winter  the  F2  region  scale  height  increases  with  magnetic  activity  (as  indicated  by  the  index). 

Apparently  this  scale  height  increase  is  the  result  of  increased  electron  and  ion  temperatures. 

It  is  not  clear  whether  heating  by  hydromagnetic  waves  (as  suggested  by  Dessler^),  by  electric 
77 

fields,  or  by  particle  precipitation  is  the  responsible  agency. 

Heating  of  the  neutral  constituent  by  an  agency  correlated  with  K  or  a  has  been  reported 
78-83  P  P 

by  several  workers.  The  time  resolution  in  these  measurements,  which  are  based  on  the 

observed  drag  on  artificial  earth  satellites,  is  inadequate  to  show  diurnal  changes.  Jacchia 
84 

and  Slowey  find  AT  /AK  ~  35 °K  for  K  ^5.  For  K  >  5,  there  is  better  correlation  with  the 
J  n7  p  p  ^  P 

a  index  and  AT  /Aa  =  1°K.  Since  our  September  results  were  obtained  at  times  when  K  was 
p  n7  p  ^  p 

mostly  <5,  we  would  expect  closest  agreement  between  these  satellite  results  and  the  values  for 
ATj/AKp  at  300  km.  The  large  difference  in  daytime  and  nighttime  values  for  this  quantity  sug¬ 
gests  that  either  (a)  the  heat  flux  associated  with  disturbed  conditions  is  absent  during  the  day¬ 
time,  or  (b)  it  cannot  be  detected  during  daytime  in  the  presence  of  the  much  larger  ultraviolet 
flux  from  the  sun.  In  either  case,  it  may  be  that  the  variation  of  T^  with  magnetic  conditions 
observed  from  satellite  results  is  largely  a  nighttime  phenomenon.  The  nature  of  this  heat 

source  remains  to  be  determined,  but  it  seems  evident  that  measurements  of  T  and  T.  over  a 

e  l 

wide  range  of  height,  made  with  a  time  resolution  far  superior  to  the  present  determinations  of 
T^,  should  contribute  significantly. 


XV.  SCALE  HEIGHT  OF  UPPER  F-REGION 


Throughout  most  of  the  daytime  hours,  the  electron-density  distribution  above  h  ex- 
hibits  almost  constant  scale  height  H  up  to  about  500  to  600  km.  The  value  of  the  scale  height 
has  been  obtained  from  the  mean  profiles  by  determining  at  what  height  the  electron  density 
falls  by  e  from  its  value  at  350  km.  This  quantity  has  a  value  of  the  order  of  100  to  150  km, 
and  is  therefore  approximately  equivalent  to  the  scale  height  at  400-  to  425 -km  altitude.  The 
values  for  the  mean  scale  height  have  been  plotted  in  Figs.  27(a)  through  (j).  No  values  are 
given  for  May  and  June  because  no  electron-density  profiles  were  obtained  in  these  months. 
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TABLE  IV 

FRACTIONAL  INCREASES  IN  TEMPERATURE  AT  VARIOUS  HEIGHTS 

Daytime 

Index 

300  km 

400  km 

500  km 

600  km 

Mean  of 

All  Heights 

AT./AA 

1  P 

1.4 

0.7 

-4.5 

-2.1 

-1.1 

AT./Aa 
'  P 

1.  1 

-0.2 

-5.1 

-2.3 

-2.2 

AT./AK 

1  P 

20 

10 

-31 

-11 

-3 

AT  /AA 
e  p 

-15.2 

-24.8 

-21.7 

-19. 1 

-20.2 

AT  /Aa 
e  p 

-12.5 

-22.0 

-19.6 

-19.8 

-18.5 

AT  /AK 
e  p 

-46 

-123 

-76 

-80 

-81 

Nighttime 

AT./AA 
'  P 

9.5 

7.6 

6.4 

- 

7.8 

AT./Aa 

1  P 

23.1 

20.9 

20.0 

- 

21.3 

AT./AK 
•  P 

195 

186 

185 

- 

189 

AT  /AA 
e  p 

15.2 

16.3 

13.2 

- 

14.9 

AT  /Aa 
e  p 

43.6 

43.  1 

45.6 

- 

44.1 

AT  /AK 
®  P 

378 

367 

422 

- 

389 
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Under  conditions  of  diffusive  equilibrium,  when  the  electron  and  ion  temperatures  and 
are  changing  with  height,  H  is  given  in 


H 


A  dogN)  =  - 


k(T.  +  Te) 


1 

T.  +  T 
1  e 


-jMT.  +  T  ) 

dz  1  e 


(5) 


in  which  z  is  geopotential  altitude  [z  =  hRe/(Re  +  h),  where  h  is  the  true  altitude  and  R0  the 

earth's  radius].  The  difference  between  true  height  h  and  geopotential  altitude  z  can  be  allowed 

_2 

for  approximately  by  taking  the  value  for  g  (868  cm  sec  at  400  km). 

+  -23 

Using  the  above  expression  and  assuming  that  O  is  the  principal  ion  (im  =  2.65  X  10  gm), 
we  have  computed  the  expected  scale  height  at  400  km  for  the  average  temperatures,  and  these 
values  are  also  shown  in  Figs.  27(a)  through  (j).  In  several  months  (e.g.,  November)  the  agree¬ 
ment  seems  excellent;  there  is  a  tendency,  however,  for  the  daytime  observed  values  to  lie 
below  the  computed  ones,  particularly  in  the  early  part  of  the  year.  It  is  possible  that  this 
arises  as  a  consequence  of  the  fact  that  the  mean  ion  mass  is  less  than  that  of  atomic  oxygen. 
However,  the  agreement  in  October  and  November  suggests  that  this  is  not  the  case,  and  a 
more  likely  explanation  is  that  in  some  months  the  density  profiles  are  in  error  as  a  result  of 
the  TR  recovery  problem  discussed  previously  (Sec.  VI -A).  In  particular,  it  is  known  that  this 
problem  was  severe  during  some  days  in  July.  It  would  seem,  therefore,  that  the  upper  part 
of  the  F-region  is  in  diffusive  equilibrium  and  that  0+  is  the  principal  ion  up  to  about  500  km 
both  by  day  and  by  night. 

This  is  an  important  result;  it  implies,  for  instance,  that  the  evening  increase  in  fQF2 
observed  in  the  summer  months  is  simply  a  consequence  of  cooling  of  the  ionosphere,  causing 
a  lowering  of  the  equilibrium  scale  height.  It  seems  quite  unnecessary  to  invoke  electrodynamic 
drift  motions  to  explain  this  phenomenon  though,  on  the  basis  of  the  results  presented  here,  one 
cannot  exclude  the  possibility  of  such  drift  motions  provided  they  are  small. 


XVI.  HEAT  FLUXQ350 

A  second  quantity  extracted  from  the  data  and  presented  in  Figs.  28(a)  through  (j)  is  the 
heat  flux  Q  defined  in 


T  -  T. 
e  l 


2  X  106QT  3/2 
e 


N 


(6) 


-3  - 1 

where  Q  is  the  energy  input  to  the  electron  gas  (ev  cm  sec  ),  and  N  is  the  electron  density. 

85  63 

This  expression  has  been  derived  by  Dalgarno,  et  al.,  and  Hanson,  and  specifies  the  heat 

input  required  to  maintain  a  given  temperature  difference  T^  —  T^  at  an  altitude  (above  300  km) 

where  the  electron  gas  is  cooled  entirely  by  coulomb  encounters  with  atomic  oxygen  ions.  Since 

T  ,  T.,  and  N  are  known  we  can  compute  Q  as  a  function  of  time  and  height.  However,  it  has 

been  shown^  that  during  the  day,  when  T  >  T.  at  all  heights,  the  altitude  variation  of  Q  is 

2  e  1 

determined  chiefly  by  the  variation  of  N  in  Eq.  (6).  That  is,  Q  decreases  with  a  scale  height 

about  half  that  for  the  electron  density.  At  night  this  is  not  necessarily  true,  and  Q  may  show 

a  peak  where  Tq/T^  reaches  a  maximum  (~350km).  Accordingly,  Q  has  been  computed  for  a 

fixed  height  of  350  km  both  for  daytime  and  nighttime.  In  the  summer  months,  Q350  is  found  to 

increase  almost  linearly  throughout  the  day  and  to  reach  a  peak  of  the  order  of  500  to  600  ev  cm" 
-1 

sec  near  sunset.  This  is  accounted  for  by  the  rise  of  the  whole  layer  throughout  the  day. 
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causing  the  electron  density  N  at  350  km  to  increase.  In  the  winter  months  the  peak  shifts  to 
the  middle  of  the  day,  as  the  density  at  350  km  is  greatest  at  that  time. 

It  should  be  stressed  that  these  daytime  values  of  Q  do  not  represent  the  total  heat  flux  into 
any  cell  of  the  electron  gas  at  350-km  altitude.  They  represent  the  heat  lost  in  such  a  cell  to 
the  ions.  A  considerably  larger  amount  of  heat  is  conducted  to  other  altitudes. 

The  flux  Q35Q  is  most  accurately  measured  when  T^  —  T^  is  large;  hence,  the  greatest  un¬ 
certainty  lies  in  values  obtained  during  the  summer  at  night.  Despite  this,  small  fluxes  in  the 
-3  -1 

range  10  to  50  ev  cm  sec  are  found  at  night  in  all  months;  there  seems  to  be  no  significant 

seasonal  variation.  These  fluxes  are  of  the  order  of  one  tenth  or  less  of  the  daytime  fluxes  but 

are  entirely  adequate  to  cause  pronounced  effects,  especially  in  winter,  when  there  is  a  large 
2 

diurnal  change  in  N  . 

It  is  possible  that  more  than  one  agency  contributes  to  this  nighttime  heating,  but  for  want 

of  better  evidence  we  presume  there  is  a  single  source  whose  intensity  is  correlated  with  K^. 

During  the  daytime,  the  peak  in  T^/T^  is  about  100  to  150  km  higher  than  the  peak  of  the  layer 

h  .  At  night,  during  those  months  where  the  T  /T.  plot  shows  a  peak,  the  height  difference 
max  e  1 

is  usually  less  than  50  km.  Evidently  the  nighttime  heating  must  occur  in  the  vicinity  of  300  km 
in  order  that  T^/T^  and  N  can  both  have  peaks  near  this  altitude.  By  contrast,  the  daytime 
absorption  of  solar  ultraviolet  takes  place  chiefly  at  heights  below  250  km. 

XVII.  DISCUSSION 

A.  High  Values  for  Electron  Temperature 

The  possible  difference  between  the  temperature  of  the  electrons  T^  and  the  neutral  gas 

temperature  T  seems  to  have  been  considered  quantitatively  only  recently.  Hanson  and 
86  ^ 

Johnson  outlined  the  physical  processes  involved,  and  their  work  has  been  summarized  by 
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Bourdeau  and  Bauer.  More  recently,  Hanson,  Dalgarno,  et  al.,  and  Dalgarno  have  dis¬ 
cussed  the  problem  anew. 

Photoelectrons  with  a  wide  range  of  energies  (5  to  50  ev)  are  produced  by  the  incident  solar 
flux.  The  first  step  in  any  calculation  of  T^  —  T^  consists  of  estimating  the  number  of  photo¬ 
electrons  produced  as  a  function  of  height  and  their  average  energy.  In  view  of  the  many  pos¬ 
sible  excitation  reactions,  this  is  perhaps  one  of  the  most  uncertain  steps  in  the  computation. 
Next,  it  is  necessary  to  investigate  the  mechanisms  by  which  the  fast  electrons  may  lose  their 

energy;  these  too  are  a  function  of  height.  At  low  levels  the  electrons  rapidly  lose  energy  by 

1 

inelastic  collisions  with  the  neutral  gas.  Important  are  the  excitation  of  the  D  state  of  atomic 

oxygen  and  also  the  excitation  of  molecular  nitrogen  into  vibrational  states.  These  effects 

should  rapidly  reduce  the  electron  energy  to  a  little  less  than  2  ev.  At  this  energy  level  a  fast 

photoelectron  is  finally  slowed  to  thermal  speeds  by  coulomb  interaction  with  other  electrons. 

This  latter  effect  raises  the  temperature  of  the  electron  gas  as  a  whole. 

Proceeding  upward  in  height,  the  number  density  of  neutral  particles  decreases  and  a  larger 

fraction  of  the  excess  energy  is  given  to  the  electron  gas  via  coulomb  encounters.  In  this  region 

(above  300  km)  the  temperature  difference  no  longer  depends  linearly  upon  the  energy  of  the 

photoelectrons,  as  those  with  more  than  about  10  ev  can  escape  along  the  field  lines  to  the  con- 
6  3 

jugate  point.  Thus,  heating  efficiency  (of  the  electron  gas  by  the  fast  electrons)  appears  to 
increase  with  height  in  a  rather  uncertain  manner.  At  these  heights  the  electron  gas  loses 
energy  primarily  to  the  (oxygen)  ions  via  elastic  collisions,  and  the  temperature  difference  will 
be  as  shown  in  Eq.  (6).^’  ^ 
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Dalgarno,  et  al.,  ignore  transport  of  heat  either  by  thermal  conduction  or  by  the  fast 

photoelectrons,  and  assume  that  the  heat  Q  given  to  the  electrons  at  any  height  is  determined 

63 

by  absorption  of  solar  ultraviolet  radiation  at  that  height.  Hanson  considers  upward  transport 

of  energetic  electrons  (which  could  contribute  considerably  to  the  value  of  Q  above  300  km)  but 

63 

is  unable  to  reach  very  firm  conclusions.  Two  other  points  which  Hanson  raises  are  (1)  the 
electrons  have  high  thermal  conductivity  and  might  be  expected  to  maintain  a  constant  tempera¬ 
ture  more  or  less  independent  of  the  atmospheric  temperature  above  600  km,  and  (2)  the  ions 
will  be  in  good  thermal  contact  with  the  neutral  particles  at  low  altitudes  but,  as  the  collision 
frequency  falls  with  height,  they  may  assume  the  electron  temperature  (above  about  750  km). 

Dalgarno^2  supports  this  latter  prediction. 

85  63 

Both  Dalgarno,  et  al.,  and  Hanson  anticipate  that  the  largest  difference  between  the 

13 

electron  and  neutral  temperatures  should  occur  near  220  km,  and  Brace,  et  al.,  have  provided 
some  experimental  evidence  in  support  of  this.  On  the  other  hand,  the  results  presented  here 
indicate  the  largest  difference  lies  between  350  and  600  km. 

In  a  more  recent  attack  on  the  problem,^  attempts  were  made  to  include  thermal  conduction 
and  an  upward  flux  of  fast  photoelectrons.  It  was  assumed  that  the  thermal  conductivity  of  the 
electrons  was  so  much  greater  than  that  of  the  ions,  that  the  latter  could  be  neglected.  Thus, 
electrons  at  heights  above  300  km  are  heated  by  (1)  heat  conducted  from  below,  and  (2)  an  up¬ 
ward  flux  of  fast  photoelectrons.  They  lose  heat  partly  by  conduction  to  other  heights  and  partly 
in  coulomb  counters  with  ions.  The  ions  are  presumed  to  lose  heat  solely  by  collisions  with  the 

neutral  particles.  This  model  has  been  able  to  account  satisfactorily  for  the  behavior  of  T  and 
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T  observed  during  the  daytime  at  midlatitudes  in  summer  at  sunspot  minimum.  Figure  29 
shows  a  comparison  of  the  temperatures  given  in  Fig.  11  with  temperatures  deduced  by  Geisler 
and  Bowhill^0  according  to  the  model.  The  agreement  is  considered  quite  good  in  view  of  the 
many  uncertain  quantities  remaining  in  the  theory.  Geisler  and  Bowhill^  report,  however, 
that  the  thermal  conductivity  of  the  electrons  is  so  great  that  many  other  parameters  (e.g.,  the 
electron-density  distribution)  can  be  changed  without  significantly  influencing  the  results.  It 
seems  therefore  that  a  theoretical  understanding  of  the  behavior  of  the  electron  and  ion  tempera¬ 
tures  in  the  ionosphere  is  now  at  hand.  This  seems  to  be  a  demonstration  of  the  power  of  the 
theoretical  method  to  provide  the  correct  answer  once  it  is  known. 

B.  Ionospheric  Anomalies 

48 

The  summer  evening  increases  in  f  F2  have  been  explained  as  a  result  of  this  work.  It  has 

been  shown  to  be  simply  a  consequence  of  the  cooling  of  the  exospheric  electron  temperature  at 
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sunset.  We  have  suggested  that  the  nocturnal  increase  in  fQF2  observed  in  winter  is  the  con¬ 
sequence  of  the  same  phenomenon  to  an  observer  in  the  winter  hemisphere.  This  latter  proposal 
lacks  a  theoretical  treatment  of  the  problem  of  temperature  distribution  along  a  field  line  when 
one  foot  is  sunlit  and  the  other  is  in  shadow.  Hence,  it  cannot  be  regarded  as  other  than  specu¬ 
lative  at  the  present  time. 

The  major  temperate  latitude  ionospheric  anomaly,  the  so-called  "winter  anomaly"  of 
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higher  midday  winter  fQF2  values  than  summer,  has  not  been  explained.  Wright  has  argued 
that  "winter  anomaly"  is  an  improper  description  since  it  can  be  shown  that  summer  behavior 
is  anomalous.  In  part,  the  seasonal  anomaly  can  be  explained  as  a  result  of  an  expansion  of  the 
F-region  in  summer,  leading  to  a  scale  height  H  above  hmax  that  is  higher  in  summer  than  in 


25 


winter.  Observations  of  Faraday  rotation  of  moon-reflected  radar  signals  yield  values  for 

go  71  72 

the  total  electron  content  of  the  ionosphere.  '  '  These  show  that  the  anomaly  in  the  total 

content  is  much  less  pronounced  than  in  the  value  for  the  peak  of  the  layer.  Nevertheless,  an 
anomaly  does  remain. 

The  results  presented  here  show  that  the  seasonal  anomaly  cannot  be  satisfactorily  ex¬ 
plained  as  a  consequence  of  large  seasonal  changes  in  the  ionospheric  temperature  which  might 
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either  (1)  violently  expand  the  region  as  suggested  by  Appleton,  or  (2)  change  the  rates  of  re- 
91 

combination.  Presumably  the  reason  for  the  anomaly  must  be  sought  in  seasonal  changes  in 

the  loss  rates  governing  the  disappearance  of  electrons  at  F-region  heights.  These  changes 
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could  be  realized  by  changing  the  composition  of  the  atmosphere.  ’ 

xvm.  CONCLUSIONS 

The  operation  of  the  Millstone  Hill  ionospheric  radar  has  been  described,  and  results  of 
about  1500  hours  of  observation  throughout  1963  have  been  presented.  Several  hitherto  un¬ 
explained  features  of  the  behavior  of  the  F-region  are  now  understood  as  a  result  of  this  work. 
The  results  reported  herein  apply  to  conditions  at  sunspot  minimum.  It  is  currently  thought 
that  conditions  should  be  substantially  different  at  sunspot  maximum.  In  particular,  it  is  be¬ 
lieved  that  the  temperature  difference  T^  —  T^  will  be  much  less.^'  ^  This  prediction  re¬ 
mains  to  be  tested.  Also,  at  sunspot  maximum  it  seems  likely  that  heating  associated  with 
magnetically  disturbed  conditions  will  be  much  more  common.  The  nature  of  the  nighttime  heat 
source  remains  to  be  established.  We  do  not  yet  know  whether  the  high  temperatures  observed 
in  September  were  the  result  of  an  increase  in  the  intensity  of  the  nighttime  agent,  or  of  yet  a 
third  source  of  ionospheric  heating. 

Many  other  opportunities  presented  by  the  development  of  the  backscatter  technique  remain 
to  be  explored.  A  more  careful  examination  of  the  uppermost  spectra  shall  be  conducted  to 
determine  the  ratio  of  the  numbers  of  0+  to  He+  ions.  This,  in  turn,  will  lead  to  temperature 
values  at  heights  which  are  less  uncertain  than  those  presented  here.  The  region  below  240  km 
remains  to  be  explored  effectively.  In  this  region,  N,  T  ,  T.  and  the  ratios  of  the  numbers  of 
ions  all  change  rapidly  with  height.  It  is  unlikely  that,  even  were  there  no  problem  with  ground 
clutter  echoes,  our  existing  radar  could  satisfactorily  examine  this  region.  (This  is  also  true 
of  the  Arecibo  Ionospheric  Observatory.)  We  intend,  therefore,  to  use  a  radar  operating  at 
1295 Mcps  with  a  beam  directed  obliquely  (at  20°  elevation)  to  explore  the  E  and  FI  regions. 

Other  groups  employing  the  incoherent  scatter  technique  exist  in  several  countries:  e.g., 
at  Prince  Albert  in  Canada;  at  Nancay  in  France;  at  Malvern  in  England;  as  well  as  the  Arecibo 
Ionospheric  Observatory  in  Puerto  Rico,  and  the  Jicamarca  Ionospheric  Radar  in  Peru.  Many 
useful  results  should  emerge  from  this  work. 
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Fig.  1.  Theoretical  spectra  of  ionospheric  backscatter  signals  computed  for  various 
values  of  electron-to-ion  temperature  ratio  Spectra  are  symmetrical  about 

center  frequency.  It  has  been  assumed  that  the  mass  m\  of  the  ions  is  that  of  0*. 
Abscissa  is  Doppler  shift  f  normalized  by  multiplying  by  radio  wavelength  X  and  a 
term  inversely  proportional  to  velocity  of  sound  for  ions. 
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Fig.  2.  Effect  on  spectrum  in  which  Te/Tj  =  2.  0  as  progressively  larger  amounts  of 
He*  ions  are  introduced  into  an  0*/He*  mixture.  Abscissa  same  as  Fig.  1 ,  and  mj 
represents  the  mass  of  the  heavier  constituent  (0*). 
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Fig.  3.  The  70-meter  parabola  employed  at  Millstone  Hill  Radar  Observatory 
for  ionospheric  studies.  Tripod  structure  supports  feedhorn  and  turnstile  junc¬ 
tion.  Left-hand  waveguide  connects  turnstile  junction  to  transmitter  (in  build¬ 
ing  off  the  picture),  and  right-hand  waveguide  connects  junction  to  firststages 
of  receiver  housed  in  lower  hut. 
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Fig.  4.  A  cathode-ray- tube  display  presenting  result  of  integration  over  a  time 
base  performed  in  computer.  In  this  measurement,  pulse  length  was  0.5  msec 
and  points  representing  averaged  output  power  are  at  100-psec  intervals.  Fol¬ 
lowing  receiver  suppression  (extreme  left)  are  first  some  ground  clutter  echoes 
and  second  the  echo  from  the  ionosphere.  Beyond  ionospheric  echo  are  echoes 
from  two  satellites  that  chanced  to  traverse  the  beam  during  this  measurement. 
Large  pulse  at  right  is  introduced  into  receiver  to  calibrate  scale  and  corre¬ 
sponds  to  an  increase  in  system  temperature  of  100°K. 
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Fig.  5.  Plot  produced  by  computer  of  electron  density  N'(h)  vs  height. 
Computer  draws  straight  lines  between  measured  points  which  are  at 
15-km  height  intervals.  Fluctuations  observed  above  500-km  height  are 
errors  due  to  noise  and  not  true  variations  in  electron  density.  Below 
200  km,  receiver  was  suppressed  against  ground  clutter  echoes. 


Fig.  6.  Example  of  combined  electron-density  profile  obtained  from 
plots  (e.  g..  Fig.  5)  for  0. 1-,  0.5-  and  1.0-msec  pulses.  These  partic¬ 
ular  measurements  were  made  at  same  time  as  rocket  was  launched  from 
Wallops  Island,  Virginia.  Difference  between  rocket  and  backscatter 
profile  is  largely  due  to  fact  that  a  has  been  assumed  constant  to  ob¬ 
tain  this  profile. 
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Fig.  7.  Variation  of  electron-to-ion  temperature  ratio  Te/T;  obtained 
from  spectra  measured  during  period  0834  to  1005  EST  on  2  July  1963. 


Fig.  8.  Backscatter  electron-density  profile  (Fig.  6)  after  correcting  for 
variation  of  electron  cross  section  a,  with  height  according  to  Eq.  (4) 
and  results  for  Te/Ti  presented  in  Fig.  7.  Rocket  measurements  were 
made  at  0922,  but  ground  clutter  echoes  obscured  F- region  peak  in 
backscatter  measurement  until  1030. 


34 


Fig.  9.  Plots  showing  variation  of  contours  of  constant  Te/Ti  (=  R)  and  1\ 
as  functionsof  scaled  quantities  f  and  x.  Ordinate  x  is  the  ratio  between 
peak  echo  power  in  the  wing  (Fig.  1)  to  that  at  center  frequency,  and  f 
is  Dopplershift  of  a  pointof  half  peak  intensity.  Dotted  lines  show  values 
in  absence  of  equipmental  effects,  whereas  solid  lines  show  change  in  x 
and  f  when  spectral  distortion  introduced  by  transmitter  pulses  and  receiver 
filters  is  allowed  for. 


f  ( kcps) 


Fig.  10.  Two  spectrum  analysis  charts  similar  to  full  lines  in  Fig.  9  are 
here  superimposed  to  illustrate  dependence  of  spectrum  shapes  upon  elec¬ 
tron  density  N  (specified  in  terms  of  a  plasma  frequency  fp).  Density  N 
enters  because  equations  contain  ratio  of  radar  wavelength  X  to  4irXp 
where  Xp  is  Debye  length  (=  ykl^/4vNe^) . 
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Fig.  11.  Variation  of  electron  (Te)  and  ion  (Tj)  temperatures  observed  on  2  July  1963. 
Solid  lines  denote  values  deduced  on  assumption  that  only  0+  ions  are  present  atall 
heights.  Dotted  lines  indicate  how  temperature  might  vary  if  percentage  of  He+  ions 
becomes  significant  above  500  km  and  increases  to  valueof  20  percentof  total  at  720  km. 


Fig.  12.  Mean  of  all  experimental  profiles  obtained  at  greatest  height  (720  km)  during 
period  0834  to  1005  on  2  July  1963  is  here  compared  with  closest  fitting  precomputed 
curves  for  an  0+  gas  and  mixture  containing  80  percent  0+,  20  percent  He+.  A  better 
fit  could  be  obtained  by  adjusting  theoretical  curves  so  that  for  0+,  Tj  =  2040°  and 
TeAi  =  1.14,  and  for  0+/He+,  Tj  =  1400°  and  Te/Ti  =  1.40.  When  this  is  done,  so¬ 
lutions  are  experimentally  indistinguishable. 
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Fig.  13.  Spectrum  analysis  chart  for  fp  =  1 .5  Mcps  and  different  ion  compositions. 
Solid  lines  represent  wholly  0+ions  in  absence  of  equipmental  effects;  dotted  lines 
denote  80/20  percent  mixture  of  0+/Fle+  ions  observed  using  a  1-msec  pulse.  If  x 
and  f  are  the  only  quantities  scaled  from  the  records,  widely  differing  values  of 
Te/Tj  can  be  obtained  depending  upon  whether  or  not  He+ ions  are  assumed  present. 


ELECTRONS /cm3 

Fig.  14.  Comparison  between  electron-density  profile  obtained  by  the 
Alouette  satellite  and  backscatter  radar  on  the  morning  of  12  July. 
Large  discrepancy  in  two  profiles  above  500  km  is  attributed  to  improper 
recovery  of  receiver  TR  at  this  time.  (The  Alouette  profile  was  kindly 
provided  by  Dr.  T.E.  Van  Zandt.) 
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Figs.  15(a-j).  Contours  of  constant  plasma  frequency  (Mcps). 
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Fig.  15.  Continued. 
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Fig.  15.  Continued. 
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Figs.  16(a-l).  Variations 
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T.  with  time  and  height. 
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Fig.  16.  Continued. 
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Fig.  16.  Continued. 
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Fig.  17.  Continued. 
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Fig.  17.  Continued. 
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Figs.  18(a-l).  Variations  of  ratio  T  /T.  with  time  and  height. 
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Fig.  18.  Continued. 
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Fig.  18.  Continued. 
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Figs.  19(a-l).  Average  daytime  behavior  of  T. ,  T^,  and  T  /T.  with  height. 
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Fig.  19.  Continued. 


51 


HEIGHT  (km)  HEIGHT  (km) 


(e) 


600  1000  1400  1800  2200  2600  3000 


TEMPERATURE  CIO 

(0 


Fig.  19.  Continued. 


52 


600  1000  1400  1800  2200  2600  3000 


TEMPERATURE  CR) 

(g) 


(h) 


Fig.  19.  Continued. 
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Fig.  19.  Continued. 
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Fig.  19.  Continued. 
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Fig.  20.  Seasonal  variations  in  average  daytime  values  of  Tj.  Solar  radio 
flux  at  10.7  cm  is  also  shown. 


Fig.  21.  Values  of  planetary  magnetic  index  Ap  on  days  when  observations  were 
conducted  in  1963,  together  with  monthly  means  of  these  values. 
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Fig.  22.  Seasonal  variations  in  daytime  Te  and  Te/Tj. 
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Figs.  23(a-l).  Average  nighttime  behavior  of  T.,  T 


and  T^/T.  with  height. 
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Fig.  23.  Continued. 
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Fig.  23.  Continued. 
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Fig.  23.  Continued. 
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Fig.  23.  Continued. 
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Fig.  23.  Continued. 
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Fig.  24.  Seasonal  variations  in  nighttime  T. ,  T  ,  and  T  /T.  at  350  km. 

i  e  e  i 


6-4 


Fig.  25.  Comparison  of  average  height  dependence  of  Te  and  Tj  on  magnetically 
quiet  days  (Ap>$  i4)  and  disturbed  days  (Ap  ^-19)  in  September.  As  can  be  seen, 
there  are  no  major  differences. 


Fig.  26.  Comparison  of  average  height  dependence  of  Te  and  T;  on  magnetically 
quiet  nights  (Ap^  14)  and  disturbed  nights  (Ap^l9)  in  September.  At  night, 
large  increases  associated  with  disturbed  conditions  occur  in  Te  and  T;. 
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Figs.  27(a-j).  Comparison  of  observed  scale  height  of  electron-density+distribution  at  400-km  altitude 
with  that  computed  from  measured  values  of  Te  and  T*,  assuming  any  O  ions  are  present. 
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Fig.  27.  Continued. 
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Fig.  27.  Continued. 
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Figs.  28(a-j).  Heat  flux  Qg5Q  required  to  maintain  temperature  difference  —  T.  observed  at  350  km. 
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Fig.  28.  Continued. 
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Fig.  28.  Continued. 


TEMPERATURE  CK ) 

Fig.  29.  Comparison  of  temperature  distribution  with  height  shown  in  Fig.  1 1 , 
with  theoretical  calculations  by  Geisler  and  Bowhill^  for  midday  in  summer 
at  sunspot  minimum. 
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